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Abstract. In this paper we construct a parametrix for the high-energy asymp- 
totics of the analytic continuation of the resolvent on a Riemannian manifold 
which is a small perturbation of the Poincare metric on hyperbolic space. As 
a result, we obtain non-trapping high energy estimates for this analytic con- 
tinuation. 



Introduction 

Under appropriate conditions, the resolvent of the Laplacian on an asymptot- 
ically hyperbolic space continues analytically through the spectrum [16] . In this 
paper we obtain estimates on the analytic continuation of the resolvent for the 
Laplacian of a metric that is a small perturbation of the Poincare metric on hy- 
perbolic space. In particular we show for these perturbations of the metric, and 
allowing in addition a real-valued potential, that there are only a finite number 
of poles for the analytic continuation of the resolvent to any half-plane containing 
the physical region and that the resolvent satisfies polynomial bounds on appropri- 
ate weighted Sobolev spaces near infinity in such a strip. This result, for a small 
strip, is then applied to the twisted Laplacian which is the stationary part of the 
d'Alembertian on de Sitter-Schwarzschild. In a companion paper [19] the decay of 
solutions to the wave equation on de Sitter-Schwarzschild space is analyzed using 
these estimates. 

In the main part of this paper constructive semiclassical methods are used to 
analyze the resolvent of the Laplacian, and potential perturbations of it, for a 
complete, asymptotically hyperbolic, metric on the interior of the ball, B" +1 , 

98 = g Q + xs(z)H, 

a) f (i-\* \y 
go = (i-izpr X5(z)=x[— r 

Here go is the standard hyperbolic metric, H = H(z, dz) is a symmetric 2-tensor 
which is smooth up to the boundary of the ball and \ G C°°(M), has x( s ) = 1 
if \s\ < |, x( s ) = if \s\ > 1. The perturbation here is always the same at the 
boundary but is cut off closer and closer to it as S 4- 0. For 6 > small enough we 
show that the analytic continuation of the resolvent of the Laplacian is smooth, so 
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has no poles, in the intersection of the exterior of a sufficiently large ball with any 
strip around the real axis in the non-physical half-space as an operator between 
weighted L 2 or Sobolev spaces, and obtain high-energy estimates for this resolvent 
in this strip. 

A special case of these estimates is as follows: Let x — , W G C°°(B n+1 ) 
be real-valued and let Rs(cr) = {A gs + x 2 W — a 2 — n/2 2 )~ 1 denote the resolvent of 
A gs + x 2 W. The spectral theorem shows that Rs(a) is well defined as a bounded 
operator in L 2 (B n+1 ; dg) if line << 0, and the results of Mazzeo and the first author 
show that it continues meromorphically to the upper half plane. Here we show that 
there exists a strip about the real axis such that if S is small and a and b are suitably 
chosen x a Rs(cr)x h has no poles, provided \a\ is large, and moreover we obtain a 
polynomial bound for the norm of x a Rs(a)x b . More precisely, with H k (B™ +1 ) the 
L 2 -based Sobolev space of order k, so for k = 0, H$(M n+1 ) = L 2 (M n+1 ;dg) : 



Theorem 1. (See Theorem M . l\ for the full statement.) There exist do > 0, such that 
*/0 < S < Sq, then x a Rs(cr)x b continues holomorphically to the region Imcr < M, 
M > \o~\ > K(5,M), provided Imcr < 6, and a > Imc. Moreover, there exists 
C > such that 

\\x a R s (a)x h v\\ HHmn+l) <CV|- 1+ *+*||i;|| £a( B»+i), fc = 0,1,2, 



\x a R s (a)x b v\\ L 2 (Kn+1) < C\a\- 1+ % +k \\v\\ H - k(Kn+1) , k = 0,1,2, 



This estimate is not optimal; the optimal bound is expected to be 0(\a\~ 1+k ). 
The additional factor of \o~\ ? results from ignoring the oscillatory behavior of the 
Schwartz kernel of the resolvent; a stationary phase type argument (which is made 
delicate by the intersecting Lagrangians described later) should give an improved 
result. However, for our application, which we now describe, the polynomial loss 
suffices, as there are similar losses from the trapped geodesies in compact sets. 

As noted above, these results and the underlying estimates can be applied to 
study the wave equation on 1 + (n + 1) dimensional de Sitter-Schwarzschild space. 
This model is given by 

(3) M = R t xX, X = [rb H ,r s i] r x§£, 
with the Lorentzian metric 

(4) G = a 2 dt 2 - a' 2 dr 2 - r 2 duj 2 , 
where duj 2 is the standard metric on §™, 

o A 9 \ 1/2 

2m AH X ' 



^ r 3 

with A and m positive constants satisfying < 9m 2 A < 1, and rbH, fsi are the two 
positive roots of a = 0. 

The d'Alembertian associated to the metric G is 

(6) □ = a- 2 (D 2 - a 2 r- n D r {r n a 2 D r ) - a 2 r - 2 A u ). 

An important goal is then to give a precise description of the asymptotic behavior, 
in all regions, of space-time, of the solution of the wave equation, \3u = 0, with 
initial data which are is necessarily compactly supported. The results given below 
can be used to attain this goal; see the companion paper [19] . 
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Since we are only interested here in the null space of the d'Alembertian, the 
leading factor of oT 2 can be dropped. The results above can be applied to the 
corresponding stationary operator, which is a twisted Laplacian 

(7) A x = a 2 r~ n D r (a 2 r n D r ) + a 2 r - 2 A u . 

In what follows we will sometimes consider a as a boundary defining function of X. 
This amounts to a change in the C°° structure of X; we denote the new manifold 
byX,. 

The second order elliptic operator Ax in ([7]) is self-adjoint, and non-negative, 
with respect to the measure 

(8) Q = a- 2 r n drdu, 

with a given by ([5]). So, by the spectral theorem, the resolvent 

(9) R(a) = (A x - a 2 )' 1 : L 2 (X; Q) — ► L 2 (X; fi) 

is holomorphic for Im a < 0. In |21j the second author and Zworski, using methods 
of Mazzeo and the first author from [16], Sjostrand and Zworski [22] and Zworski 
[25] prove that the resolvent family has an analytic continuation. 

Theorem 2. (Sd Barreto- Zworski, see [21] ) As operators 

R(a) :C °°(i) — >C°°(i) 

the family (|9]) has a meromorphic continuation to C with isolated poles of finite 
rank. Moreover, there exists e > such that the only pole of R(a) with Imcr < e is 
at a = 0; it has multiplicity one. 

Theorem [2] was proved for n + 1 = 3, but its proof easily extends to higher dimen- 
sions. 

In order to describe the asymptotics of wave propagation precisely on M via 
R(cr) it is necessary to understand the action of R(a) on weighted Sobolev spaces 
for a in a strip about the real axis as | Re <j\ — > oo. The results of [TB] actually show 
that R(a) is bounded as a map between the weighted spaces in question; the issue 
is uniform control of the norm at high energies. The strategy is to obtain bounds 
for R(a) for Re(er) large in the interior of X, then near the ends TbH and r s i, and 
later glue those estimates. In the case n + 1 = 3, we can use the following result of 
Bony and Hafner to obtain bounds for the resolvent in the interior 

Theorem 3. (Bony-Hdfner, see [1]^ There exists e > and M > such that if 

o 

Imo~\ < e and \ Recr| > 1, then for any \ e Cfi°(X) there exists C > such that if 
n = 3 in then 

(10) \\xR(o-)xf\\LHX;n) < C\a\ M \\f\\ L2{X; ny 

This result is not known in higher dimensions (though the methods of Bony and 
Hafner would work even then), and to prove our main theorem we use the results 
of Datchev and the third author [5] and Wunsch and Zworski [24, to handle the 
general case. The advantage of the method of [5] is that one does not need to 
obtain a bound for the exact resolvent in the interior and we may work with the 
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approximate model of [24] instead. We decompose the manifold X in two parts 

X = X U Xi, where 
(11) X = [r bH , r bH + 4(5) x §" U (r s j - 45, r B i] x §" and 

X 1 = (r hH + 5 1 r sl -5)xS n . 

If <5 is small enough and if -f(t) is an integral curve of the Hamiltonian of Ax then 
(see Section [8] and either [5] or [14] ^) 

, , u „ dx(t) d 2 x(t) 
12 if x(j(t)) < 45 and — ^ = => — ^ < 0. 

We consider the operator Ax restricted to X%, and place it into the setting of 
[2"1] as follows. Let X[ be another Riemannian manifold extending X x ~ (r bH - + 
<5/2, r s i — (5/2, r s i) x S™ (and thus Xi) and which is Euclidean outside some compact 
set, and let Axj be a self-adjoint operator extending Ax with principal symbol 
given by the metric on X[ which is equal to the Euclidean Laplacian on the ends. 
Let 

(13) Pi = h 2 A x[ -iT, he (0, 1) 

where T 6 C°°(X[; [0, 1]) is such that T = on X x and T = 1 outside X x . Thus, 
Pi — 1 is semiclassically elliptic on a neighborhood of X[ \X\. (In particular, this 
implies that no bicharacteristic of Pi — 1 leaves X\ and returns later, i.e. X\ is 
bicharacteristically convex in X[, since this is true for X\ inside X\, and Pi — 1 
is elliptic outside X\, hence has no characteristic set there.) By Theorem 1 of [24] 
there exist positive constants c, C and e independent of h such that 

(14) \\{Pi-<J 2 )- 1 \\ L 2^ L 2<Ch- N 1 ae (1- c, 1 + c) x {-c,eh) C C. 

Due to the fact that < 9to 2 A < 1, the function /3(r) — \^a 2 {~r) satisfies 
/3(n3H) > and /3(r s i) < 0. Set /3 b H = /?(r b H) and (3 S \ = fi{r s \). The weight function 
we consider, a £ C°([r b H, f s i\), is positive in the interior and satisfies 



(15) a{r) 



a i/0bH near r bH - 
qA/IAi near r s \ 



We will prove 

Theorem 4. If n — 2, fe£ e > &e smc/i i/iai (|10p holds. In general, assume that 5 
is such that (|12p is satisfied, and let e > 6e smc/i £/iai (| 14|) holds. If 

< 7 < min(e,^ 6ff , |/3 S /|, 1), 

i/ien /or 6 > 7 i/iere exist C and M swc/i i/iai iflmcr < 7 and | Recr| > 1, 

(16) ||5 b P((7)(5 b /ll^(A;f2) < C|(7| M ||/|| L 2 ( x ; o), 

where Q is defined in ([8]). 

This result can be refined by allowing the power of the weight, on either side, to 
approach Imer at the expense of an additional logarithmic term, see Theorem llO.il 

Two proofs of Theorem |4] are given below. The first, which is somewhat simpler 
but valid only for n + 1 = 3, is given in sections M and EH It nses techniques of 
Bruneau and Petkov [5] to glue the resolvent estimates from Theorem [T] and the 
localized estimate (ITU)) . The second proof, valid in general dimension, uses the 
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estimate (fM)) and the semiclassical resolvent gluing techniques of Datchev and the 
third author [5] . This is carried out in section [TTJ 

Related weighted L 2 estimates for the resolvent on asymptotically hyperbolic 
and asymptotically Euclidean spaces have been proved by Cardoso and Vodev in 
[3]. However, such estimates, combined with Theorem [3] only give the holomorphic 
continuation of the resolvent, as an operator acting on weighted L 2 spaces, for 
Rec > 1, to a region below a curve which converges polynomially to the real axis. 
These weaker estimates do suffice to establish the asymptotic behavior of solutions 
of the wave equation modulo rapidly decaying terms (rather than exponentially 
decaying) and would give a different proof of the result of Dafcrmos and Rodnianski 

In the case of a non-trapping asymptotically hyperbolic manifold which has 
constant sectional curvature near the boundary, it was shown by Guillarmou in [9] 
that there exists a strip about the real axis, excluding a neighborhood of the origin, 
which is free of resonances. In the case studied here, the sectional curvature of the 
metric associated to Ax is not constant near the boundary, and there exist trapped 
trajectories. However, see [21], all trapped trajectories of the Hamilton flow of Ax 
are hyperbolic, and the projection of the trapped set onto the base space is contained 
in the sphere r = 3m, which is known as the ergo-sphere. Since the effects of the 
trapped trajectories are included in the estimates of Bony and Hafner, constructing 
the analytic continuation of the resolvent of a twisted Laplacian that has the correct 
asymptotic behavior at infinity, uniformly at high energies, allows one to obtain the 
desired estimates on weighted Sobolev spaces via pasting techniques introduced by 
Bruneau and Petkov [2] . The main technical result is thus Theorem [T] and its 
strengthening, Theorem ll.il 

Theorem 1 1.1 1 is proved by the construction of a high-energy parametrix for A gs + 
x 2 W. More precisely, as customary, the problem is translated to the construction 
of a semiclassical parametrix for 

P(h, a) = h 2 A g + h 2 x 2 W -h 2 ^-a 2 = h 2 (a 9 + x 2 W - ^ - (£) 

where now a G (1 — c, 1 + c) x {—Ch, Ch) C C, c, C > 0, and h G (0, 1), h ->• 0, so 
the actual spectral parameter is 

n 2 /ex 2 

T + (h) • 

and Im %■ is bounded. Note that for Im a < 0, 

(17) R(h,a) = PiKa)- 1 : L 2 {M n+1 ) ^ H 2 (B n+1 ) 

is meromorphic by the results of Mazzeo and the first author [16]. Moreover, while 
a is not real, its imaginary part is 0(h) in the semiclassical sense, and is thus not 
part of the semiclassical principal symbol of the operator. 

The construction proceeds on the semiclassical resolution Mo .h of the product of 
the double space B™ +1 x B™ +1 introduced in [16] . and the interval [0, 1)^ - this space 
is described in detail in Section |3] Recall that, for fixed h > 0, the results of [Tfi] 
show that the Schwartz kernel of P(h, cr)" 1 , defined for Imcr < 0, is well-behaved 
(polyhomogeneous conormal) on B n+1 XoB" +1 , and it extends meromorphically 
across Imcr = with a similarly polyhomogeneous conormal Schwartz kernel. Thus, 
the space we are considering is a very natural one. The semiclassical resolution is 
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already needed away from all boundaries; it consists of blowing up the diagonal at 
h = 0. Note that P(h, a) is a semiclassical differential operator which is elliptic 
in the usual sense, but its semiclassical principal symbol g — (Reer) 2 is not elliptic 
(here g is the dual metric function). Ignoring the boundaries for a moment, when 
a semiclassical differential operator is elliptic in both senses, it has a parametrix in 
the small semiclassical calculus, i.e. one which vanishes to infinite order at h = 
off the semiclassical front face (i.e. away from the diagonal in B™ +1 x B™ +1 x {0}). 
However, as P(h, a) is not elliptic semiclassically, semiclassical singularities (lack 
of decay as h — > 0) flow out of the semiclassical front face. 

It is useful to consider the flow in terms of Lagrangian geometry. Thus, the 
small calculus of order — oo semiclassical pseudodiffcrcntial operators consists of 
operators whose Schwartz kernels are semiclassical-conormal to the diagonal at 
h = 0. As P(h, a) is not semiclassically elliptic (but is elliptic in the usual sense, 
so it behaves as a semiclassical pseudodifferential operator of order — oo for our 
purposes), in order to construct a parametrix for P(h, er), we need to follow the 
flow out of semiclassical singularities from the conormal bundle of the diagonal. 
For P(h, a) as above, the resulting Lagrangian manifold is induced by the geodesic 
flow, and is in particular, up to a constant factor, the graph of the differential of 
the distance function on the product space. Thus, it is necessary to analyze the 
geodesic flow and the distance function; here the presence of boundaries is the main 
issue. As we show in Section [2] the geodesic flow is well-behaved on M n+1 Xq M n+1 
as a Lagrangian manifold of the appropriate cotangent bundle. Further, for S > 
small (this is where the smallness of the metric perturbation enters), its projection 
to the base is a diffeomorphism, which implies that the distance function is also 
well-behaved. This last step is based upon the precise description of the geodesic 
flow and the distance function on hyperbolic space, see Section [2j 

In Section[5j we then construct the parametrix by first solving away the diagonal 
singularity; this is the usual elliptic parametrix construction. Next, we solve away 
the small calculus error in Taylor series at the semiclassical front face, and then 
propagate the solution along the flow-out by solving transport equations. This 
is an analogue of the intersecting Lagrangian construction of the first author and 
Uhlmann [18], see also the work of Hassell and Wunsch [IT] in the semiclassical 
setting. So far in this discussion the boundaries of Mq^ arising from the boundaries 
of B n+1 Xo B n+1 have been ignored; these enter into the steps so far only in that 
it is necessary to ensure that the construction is uniform (in a strong, smooth, 
sense) up to these boundaries, which the semiclassical front face as well as the lift 
of {h = 0} meet transversally, and only in the zero front face, i.e. at the front face 
of the blow-up of B" +1 x B" +1 that created B" +1 x B n+1 . Next we need to analyze 
the asymptotics of the solutions of the transport equations at the left and right 
boundary faces of B™ +1 Xo B™ +1 ; this is facilitated by our analysis of the flowout 
Lagrangian (up to these boundary faces). At this point we obtain a parametrix 
whose error is smoothing and is 0(h°°), but does not, as yet, have any decay at the 
zero front face. The last step, which is completely analogous to the construction of 
Mazzeo and the first author, removes this error. 

As a warm-up to this analysis, in Section |4] we present a three dimensional 
version of this construction, with worse, but still sufficiently well-behaved error 
terms. This is made possible by a coincidence, namely that in K 3 the Schwartz 
kernel of the resolvent of the Laplacian at energy (A — z0) 2 is a constant multiple 
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of r , and r" 1 is a homogeneous function on R 3 , which enables one to blow- 
down the semiclassical front face at least to leading order. Thus, the first steps 
of the construction are simplified, though the really interesting parts, concerning 
the asymptotic behavior at the left and right boundaries along the Lagrangian, are 
unchanged. We encourage the reader to read this section first as it is more explicit 
and accessible than the treatment of arbitrary dimensions. 

In Section |6l we obtain weighted L 2 -bounds for the parametrix and its error. In 
Section [7] we used these to prove Theorem [T] and Theorem ll.il 

In Section [3] we describe in detail the de Sitter-Schwarzschild set-up. Then in 
Section [9j in dimension 3 + 1, we describe the approach of Bruneau and Pctkov 
[2] reducing the necessary problem to the combination of analysis on the ends, i.e. 
Theorem ll.il and of the cutoff resolvent, i.e. Theorem [3] Then, in Section [101 wc 
use this method to prove Theorem[4l and its strengthening, Theorem llO.il Finally, 
in Section [11] we give a different proof which works in general dimension, and does 
not require knowledge of estimates for the exact cutoff resolvent. Instead, it uses 
the results of Wunsch and Zworski [24] for normally hyperbolic trapping in the 
presence of Euclidean ends and of Datchev and the third author 5, which provide a 
method to combine these with our estimates on hyperbolic ends. Since this method 
is described in detail in [5], we keep this section fairly brief. 



1. Resolvent estimates for model operators 

In this section we state the full version of the main technical result, Theorem [T] 
Let go be the metric on M n+1 given by 

Adz 2 

(1-1) 90- (1 _| Z |2)2- 

We consider a one-parameter family of perturbations of go supported in a neigh- 
borhood of <9B n+1 of the form 

(1.2) g s = g + xs(z)H(z,dz), 

where H is a symmetric 2-tensor, which is C°° up to dW l+1 . \ G C°°(K), with 
X (s) = 1 if \s\ < i, X (s) = if |«| > 1, and Xs (z) = x ((l - M)^ 1 ). 

Let x = j^jf}, W G C°°(B" +1 ) and let R s (a) = {A gs + x 2 W - a 2 - l 2 )" 1 denote 

the resolvent of A gs +x 2 W. The spectral theorem gives that Rs(a) is well defined as 
a bounded operator in L 2 (B n+1 ) = L 2 (M n+1 ;dg) if Ima « 0. The results of [16] 
show that Rs{cr) continues meromorphically to C\«N/2 as an operator mapping C°° 
functions vanishing to infinite order at dM n+1 to distributions in B n+1 . We recall, 
see for example [TS], that for k G N, 

H*(M n+1 ) = {ue L 2 (B n+1 ) : (xd x , d^) m u G i 2 (B n+1 ), m < k}, 

and 

H- k = {»£ P'(B" +1 ) : there exists up £ L 2 (V> n+1 ),v = ^ (xd x , d^fup}. 

\p\<k 

Our main result is the following theorem: 

Theorem 1.1. There exist So > 0, such that if < 6 < So, then x a Rs{a)x b 
continues holomorphically to Ima < M, M > 0, provided \a\ > K(5,M), b > Ima 



8 



RICHARD MELROSE, ANTONIO SA BARRETO AND ANDRAS VASY 



and a > Im a. Moreover, there exists C > such that 

\\x a R s (a)x b v\\ H , iBn+1) <C\a\- 1+ ^ +k \\v\\ L 2 (Bn+1} , fc = 0,1,2, 

' ' \\x a R s (a)x b v\\ L 2 (Bn+1) < C\a\~ 1+ ^ +k \\v\\ H - k{Bn+1) , k = 0,1,2, 

If a = Imcr, or b — Imcr, or a = b = Imcr, Zei 4>n{x) £ C°°((0, 1)), 0jv > 1, 
4>n(x) = \logx\~ N , if x < j 4>n(x) = 1 ?/ a; > |. TTien m each case the operator 

T a ,b,N(cr) =■ x lma '4>N(x)Rs(o-)x b , i/ b > Imcr, 

(1.4) T^Ar^X^^Cr)^" 1 ^^), if 0, > IUXCT, 

T a .b,N =: ^"^(^^(cr)^'""^^), 

continues holomorphically to Imcr < M, provided N > i, |cr| > K(S,M). Moreover 
in each case there exists C — C{M,N,S) such that 

l|r a , 6 Jv(crH| Hofe(B „ +1) <CV|- 1 + f + fe |M| i2(B „ + 1) , fc = 0,l,2, 
||T , a,fc : A(cr)w||L2( B „ + i ) <C|cr| _ T+ ' 1 1«| | , = 0, 1, 2. 

2. The distance Function 

In the construction of the uniform parametrix for the resolvent we will make 
use of an appropriate resolution of the distance function, and geodesic flow, for 
the metric gs- This in turn is obtained by perturbation from 8 = 0, so we start 
with an analysis of the hyperbolic distance, for which there is an explicit formula. 
Namely, the distance function for the hyperbolic metric, go , is given in terms of the 
Euclidean metric on the ball by 

dist : (B")° x (B" )° — > M where 

, , m 2\z-z'\ 2 

cosh(disto(^, z )) = 1 + 



(l-|z| 2 )(l-|z'| 2 )' 

We are particularly interested in a uniform description as one or both of the 
points approach the boundary, i.e. infinity. The boundary behavior is resolved 
by lifting to the 'zero stretched product' as is implicit in [15]. This stretched 
product, B" XoB™, is the compact manifold with corners defined by blowing up the 
intersection of the diagonal and the corner of B" x B" : 

/3 : B n x B n = [B n x B"; <9Diag] — > 1" x B", 

(2.2) <9Diag = Diagn (3B™ x <9B") = {(z,z); \z\ = 1}, 

Diag = {(z, z') e B" x B"; z = z'}. 

See Figure [T] in which (x,y) and (x\y') are local coordinates near a point in the 
center of the blow up, with boundary defining functions x and x' in the two factors. 

Thus B™ x B™ has three boundary hypersurfaces, the front face introduced by 
the blow up and the left and right boundary faces which map back to dM n x B™ and 
B" x SB" respectively under (3. Denote by Diag the lift of the diagonal, which in 
this case is the closure of the inverse image under /3 of the interior of the diagonal 
in B™ x B™. 

Lemma 2.1. Lifted to the interior of B" x B™ the hyperbolic distance function 
extends smoothly up to the interior of the front face, in the complement o/Diag , 
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x 




Figure 1. The stretched product B™ x B™. 



where it is positive and, for an appropriate choice pl G C 00 (B™ Xo B") of defining 
function for the left boundary and with pn its reflection, 



with F 2 a quadratic defining function for Diag . 

Proof. We show first that the square of the Euclidean distance function, \z — z'\ 2 , 
lifts to be smooth on B™ Xo B™ and to vanish quadratically on Diag and on the 
front face produced by the blow up 



Here / > vanishes precisely at Diag and does so quadratically. Indeed this is 
certainly true away from the front face produced by the blow up. The spaces and 
the distance function are invariant under rotational symmetry, which lifts under 
the blow up, so me may fix the spherical position of one variable and suppose that 
z' = (1 — x' , 0), with x' > and small, the blow up is then of z — z', x' — 0. The 
fact that the variables are restricted to the unit ball is now irrelevant, and using 
the translation-invariance of the Euclidean distance we can suppose instead that 
z' = (x',0) and blow up z — z', x 1 — 0. Since \(x, 0) — z\ 2 is homogeneous in all 
variables it lifts to be the product of the square of a defining function for the front 
face and a quadratic defining function for the lift of the diagonal. This proves (|2.4p 
after restriction to the preimage of the balls and application of the symmetry. 

The hyperbolic distance is given by (I2.1[) where 1 — \ z\ 2 and l — \z'\ 2 are boundary 
defining functions on the two factors. If R is a defining function for the front face 
of B" Xq B™ then these lift to be of the form plR and prR so combining this with 



(2.3) 



/3*dist (z, z) = - \og(p L p R ) + F, 
< F G C°°(B" x B" \ Diag ), F 2 G C°°(B" x M n ) 



(2.4) 



(3*(\ z - z '\ 2 ) = R 2 f, f G C°°(B" x B n ). 



(E3) 



(2.5) 



(3* cosh(dist (z,z')) = 1 + 2 



PlPr 



Now exp(£) = cosht + [cosh 2 t — f] 2 , for t > 0, and from (|2.5[) it follows that 



(2.6) 
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Near Diag the square-root is dominated by the first part and near the left and 
right boundaries by the second part, and is otherwise positive and smooth. Taking 
logarithms gives the result as claimed, with the defining function taken to be one 
near Diag and to be everywhere smaller than a small positive multiple of (1 — 

\z?)/R- ' □ 

This result will be extended to the case of a perturbation of the hyperbolic met- 
ric by constructing the distance function directly from Hamilton- Jacobi theory, i.e. 
by integration of the Hamilton vector field of the metric function on the cotangent 
bundle. The presence of only simple logarithmic singularities in (|2.3[) shows, per- 
haps somewhat counter-intuitively, that the Lagrangian submanifold which is the 
graph of the differential of the distance should be smooth (away from the diagonal) 
in the b-cotangent bundle of M5 2 . Conversely if this is shown for the perturbed 
metric then the analogue of (|2.3[) follows except for the possibility of a logarithmic 
term at the front face. 

Since the metric is singular near the boundary, the dual metric function on 
T*M n is degenerate there. In terms of local coordinates near a boundary point, x, 
y where the boundary is locally x — 0, and dual variables £, 77, the metric function 
for hyperbolic space is of the form 

(2.7) 2p = x 2 i 2 + Ax 2 (l - x 2 )~ 2 h Q (uj, 77) 

where ho is the metric function for the induced metric on the boundary. 

Recall that the 0-cotangent bundle of a manifold with boundary M, denoted 
°T*M, is a smooth vector bundle over M which is a rescaled version of the ordinary 
cotangent bundle. In local coordinates near, but not at, the boundary these two 
bundles are identified by the (rescaling) map 

(2.8) T*M 3 (x, y, £, 77) ^ (x, y, A, /x) = (x, y, a& 2:77) g °T*M. 

It is precisely this rescaling which makes the hyperbolic metric into a non-degenerate 
fiber metric, uniformly up to the boundary, on this bundle. On the other hand the b- 
cotangent bundle, also a completely natural vector bundle, is obtained by rescaling 
only in the normal variable 

(2.9) T*M 3 (x, y, £, 7,) .— ► {x, y, A, 77) = (x, y, x£, 77) g b T*M . 
Identification over the interior gives natural smooth vector bundle maps 

(2.10) i b0 : °TM — ► h TM, i l h0 : h T*M — > °T*M. 

The second scaling map can be constructed directly in terms of blow up. 

Lemma 2.2. If T*dM C °Tg M M denotes the annihilator of the null space, over 

the boundary, of bbo in (|2.10p then there is a canonical diffeomorphism 

(2.11) 

h T*M — > [ T*M,°T*dM] \ f3* (°T£ M M) , f3 : [°T*M,°T*dM] — > °T*M, 

to the complement, in the blow up, of the lift of the boundary: 

f3*(°T% M M) = f3-i(°T* M M\°T*dM). 

Proof. In local coordinates x, yj, the null space of tbo in (|2.10[) is precisely the span 
of the 'tangential' basis elements xd Vj over each boundary point. Its annihilator, 
°T*dM is given in the coordinates by fj, = at x = 0. The lift of the 'old 

boundary' x — is precisely the boundary hypersurface near which \n\ dominates 
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x. Thus, x is a valid defining function for the boundary of the complement, on the 
right in (|2.11j) and locally in this set, above the coordinate patch in M, r)j = fJ,j/x 
are smooth functions. The natural bundle map h T*M — > °T*M underlying (|2.11[) 
is given in these coordinates by X^- + r\ ■ dy i — > A— + xn ■ ™, which is precisely 
the same map, p = xrj, as appears in ()2.9|) . so the result, including naturality, 
follows. □ 

The symplectic form lifted to °T*M is 

°uj = —dX A dx H — d/i A dy ^dx A (u ■ dy) 

X X x z 

whereas lifted to h T*M it is 

(2.12) b ui — —dX Adx + drj A dy. 

x 

Working, for simplicity of computation, in the non-compact upper half-space 
model for hyperbolic space the metric function lifts to the non-degenerate quadratic 
form on °T*M : 

(2.13) 2 Po = A 2 + /i (w,/i) 

where ho = |/x| 2 is actually the Euclidean metric. The O-Hamilton vector field 
of po E C°° (°T* M), just the lift of the Hamilton vector field over the interior, is 
determined by 

(2.14) %(-, H p ) = dp. 
Thus 

"p- x dX a *+ x dfl °v dp dx)° x \ 9A M+ dy) 

and hence 

(2.15) °H P0 = X{x8 x + ndp) - h d x + |H fco 

is tangent to the smooth (up to the boundary) compact sphere bundle given by 
Po = 1. 

Over the interior of M = M n , the hyperbolic distance from any interior point 
of the ball is determined by the graph of its differential, which is the flow out 
inside pq = 1, of the intersection of this smooth compact manifold with boundary 
with the cotangent fiber to the initial point. Observe that H po is also tangent to 
the surface /i = over the boundary, which is the invariantly defined subbundlc 
°T*dM. Since the coordinates can be chosen to be radial for any interior point, it 
follows that all the geodesies from the interior arrive at the boundary at x = 0, 
fi = 0, corresponding to the well-known fact that hyperbolic geodesies are normal to 
the boundary. This tangency implies that H po lifts under the blow up of °T*dM in 
(|2.1ip to a smooth vector field on h TM; this can also be seen by direct computation. 

Lemma 2.3. The graph of the differential of the distance function from any interior 
point, p E (B")°, of hyperbolic space extends by continuity to a smooth Lagrangian 
submanifold of b T*(M n \ {p}) which is transversal to the boundary, is a graph over 
B™ \ {p} and is given by integration of a non-vanishing vector field up to the bound- 
ary. 
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Proof. Observe the effect of blowing up [i = 0, x — on the Hamilton vector field 
in (|2.15|) . As noted above, near the front face produced by this blow up valid 
coordinates are given by jj — (J,/x, A and y, with x the boundary defining function. 
Since this transforms °uj to h uj it follows that °H po is transformed to 

(2.16) h H po = x(\d x - xh d x + xH ho ) 

where now H/j is the Hamilton vector field with respect to y and rj. 

The constant energy surface po = 1 remains smooth, but non-compact, near the 
boundary, which it intersects transversally in A = 1. From this the result follows 
- with the non- vanishing smooth vector field being b H po divided by the boundary 
defining function x near the boundary. □ 

The logarithmic behaviour of the distance function in (|2 .3[) with one point fixed 
in the interior is a consequence of Lemma 12.21 since the differential of the distance 
must be of the form adx/x + b ■ dy for smooth functions a and 6, and since it is 
closed, a is necessarily constant on the boundary. 

To examine the distance as a function of both variables a similar construction 
for the product, in this case M 2 , M = B" can be used. The graph, A, of the 
differential of the distance d(p,p') as a function on M 2 is the joint flow out of the 
conormal sphere bundle to the diagonal in T*M 2 = T*M x T*M, under the two 
Hamiltonian vector fields of the two metric functions within the product sphere 
bundles. As before it is natural to lift to °T*M x °T*M where the two sphere 
bundles extend smoothly up to the boundary. However, one can make a stronger 
statement, namely that the lifted Hamilton vector fields are smooth on the b- 
cotangent bundle of Mq, and indeed even on the 'partially b'-cotangent bundle of 
Mq, with 'partially' meaning it is the standard cotangent bundle over the interior 
of the front face. This is defined and discussed in more detail below; note that the 
identification of these bundles over the interior of M 2 extends to a smooth map 
from these bundles to the lift of °T*M x °T*M, as we show later, explaining the 
'stronger' claim. 

Smoothness of the Hamilton vector field together with transversality conditions 
shows that the flow-out of the conormal bundle of the diagonal is a smooth La- 
grangian submanifold of the cotangent bundle under consideration; closeness to a 
particular Lagrangian (such as that for hyperbolic space) restricted to which pro- 
jection to Mq is a diffeomorphism, guarantees that this Lagrangian is also a graph 
over Mq. Thus, over the interior, (Ad 2 ) , it is the graph of the differential of the 
distance function, and the latter is smooth; the same would hold globally if the La- 
grangian were smooth on T* Mq. The latter cannot happen, though the Lagrangian 
will be a graph in the b-, and indeed the partial b-, cotangent bundles over Mq. 
These give regularity of the distance function, namely smoothness up to the front 
face (directly for the partial b bundle, with a short argument if using the b bundle) , 
and the logarithmic behavior up to the other faces. Note that had we only showed 
the graph statement in the pullback of °T*M x °T*M, one would obtain directly 
only a weaker regularity statement for the distance function; roughly speaking, the 
closer the bundle in which the Lagrangian is described is to the standard cotangent 
bundle, the more regularity the distance function has. 

In fact it is possible to pass from the dual of the lifted product 0-tangent bundle 
to the dual of the b-tangent bundle, or indeed the partial b-bundle, by blow-up, as 
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for the single space above. Observe first that the natural inclusion 

(2.17) i oh x i b : °TM x °TM — ► h TM 2 = h TM x h TM 

identifies the sections of the bundle on the left with those sections of the bundle 
on the right, the tangent vector fields on M 2 , which are also tangent to the two 
fibrations, one for each boundary hypersurface 

(2.18) 4> L ■ dM x M — > M, (f> R : M x dM — > M. 

Lemma 2.4. The fibrations (|2.18[) , restricted to the interiors, extend by continuity 
to fibrations <pL, resp. 4>r, of the two 'old' boundary hypersurfaces of Mq and the 
smooth sections of the lift of °TM x °TM to M 2 are naturally identified with the 
subspace of the smooth sections of h TM$ which are tangent to these fibrations and 
also to the fibres of the front face of the blow up, /3o : S(Mq) — > dM x dM. 

Proof. It is only necessary to examine the geometry and vector fields near the front 
face produced by the the blow up of the diagonal near the boundary. Using the 
symmetry between the two factors, it suffices to consider two types of coordinate 
systems. The first is valid in the interior of the front fact and up to a general point 
in the interior of the intersection with one of the old boundary faces. The second 
is valid near a general point of the corner of the front face, which has fibers which 
are quarter spheres. 

For the first case let x, y and x' , y' be the same local coordinates in two factors. 
The coordinates 

(2.19) s = x/x', x', y and Y = (y' - y)/x' 

are valid locally in Mq above the point x = x' = 0, y — y', up to the lift of the old 
boundary x = 0, which becomes locally s = 0. The fibration of this hypersurface is 
given by the constancy of y and the front face is x' = with fibration also given 
by the constancy of y. The vector fields 

xd x , xd y , x d x i and x'd y > 

lift to 

sd s , sx'dy — sdy, x 1 d x i — sd s — Y-dy and dy. 

The basis sd s , sx'd y , x'd x > and dy shows that these vector fields are locally precisely 
the tangent vector fields also tangent to both fibrations. 

After relabeling the tangential variables as necessary, and possibly switching 
their signs, so that y[ — y\ > is a dominant variable, the coordinate system 

/ x x' y'j - yj 

(2.20) t = y l -yi, si = , s 2 = — , Zj = — , j > 1, y 

Vi -yi Vi - vi Vi - vi 

can be used at a point in the corner of the front face. The three boundary hyper- 
surfaces are locally s\ — 0, s 2 = and t — and their respective fibrations are 
given in these coordinates by 

si = 0, y = const., 

(2.21) s 2 = 0, y[ = y\ + 1 = const., y'j = y 3 + tZj = const., j > 1, 

t = 0, y = const. 

Thus, the intersections of fibres of the lifted left or right faces with the front face 
are precisely boundary hypersurfaces of fibres there. On the other hand within the 
intersection of the lifted left and right faces the respective fibres are transversal 
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except at the boundary representing the front face. The lifts of the basis of the 
zero vector fields is easily computed: 

xd x i — > sid Sl , 
xd yi — > -s\td t + s\d Sl + sis 2 d S2 + s\Z ■ d z + s\td Vl , 
xdy 3 i — > sitd yj - sid Zj , 
x'd x i i — > s 2 d s 



(2.22) 



J S2 J 

„2; 



x'dy^ — ► s 2 td t - S2S!d Sl - sld S2 - s 2 Z ■ d z , 



x'd y >. i — > s 2 d Zj ■ 

The span, over C°°(Mg), of these vector fields is also spanned by 

sid Sl , s 2 d S21 s 2 d Z] , j > 1, s 1 (td Vj -d Zj ), j > 1, s 1 (td t -td Vl -Z-d z ) and s 2 td t . 

These can be seen to locally span the vector fields tangent to all three boundaries 
and corresponding fibrations, proving the Lemma. □ 



With <f) — {4>l, 4>Ri A)} the collection of boundary fibrations, we denote by ^TMq 
the bundle whose smooth sections are exactly the smooth vector fields tangent to 
all boundary fibrations. Thus, the content of the preceding lemma is that 

(3*(°TM x °TM) = ^TMq. 

These fibrations allow the reconstruction of b T*Mg as a blow up of the lift of 
°T*M x °T*M to Mq. It is also useful, for more precise results later on, to consider 
the 'partially b-' cotangent bundle of Mq, b ' T*Mq; this is the dual space of the 
partially b-tangent bundle, ' S TM§ , whose smooth sections are smooth vector fields 
on Ml which are tangent to the old boundaries, but not necessarily to the front 
face, ff. Thus, in coordinates (|2.19[) . sd s , d x i, d y and dy form a basis of ' TMq, 
while in coordinates (I2.20[) . Sid si , s 2 d S2 , d t , d Zj and d y do so. Let 

i 2 Qh : ^TMq -4- b TM 2 , ^ b)ff : ^TA/ 2 bfl TM 2 
be the inclusion maps. 

Lemma 2.5. The annihilators, in the lift of°T*M x°T*M to Mq, of the null space 
of either Lg b or t^ b ff over the old boundaries, as in Lemma \2.4\ form transversal 
embedded p-submanifolds. After these are blown up, the closure of the annihilator 
of the nullspace of L^ b , resp. tg&ff' over the interior of the front face of Mq is a 
p-submanifold, the subsequent blow up of which produces a manifold with corners 
with three 'old' boundary hyper surf aces; the complement of these three hyper surf aces 
is canonically diffeomorphic to b T*M$, resp. b ' T*Mq. 

Proof. By virtue of Lemma f2.2l and the product structure away from the front face ff 
of M§, the statements here are trivially valid except possibly near ff . We may again 
use the coordinate systems discussed in the proof of Lemma l2~4l Consider the linear 
variables in the fibres in which a general point is lxd x + v ■ xd y + l'x'd x i + v' ■ x'd y i . 

First consider the inclusion into b TM 2 . In the interiors of s% = and s 2 = 
and the front face respectively, the null bundles of the inclusion into the tangent 
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vector fields are 
(2.23) 

1 = 1' = 0, v' = 0, 
1 = 1' = 0, v = 0, 

lsd s + v(sx'd y - sd Y ) + l'(x'd x > - sd s - Yd Y ) + v'dy = Q at x' = 0, s > 

<=>l = l' = 0, v' = sv. 

The corresponding annihilator bundles, over the interiors of the boundary hy- 
persurfaces of Mg , in the dual bundle, with basis 

, dx dy .,dx' ,dy' 

(2.24) A h fi— + A — - + M -T" 

a: x a;' a;' 

are therefore, as submanifolds, 

si = 0, /tx = 0, 

(2.25) s 2 =0, fi' = and 

a;' = 0, /i + s/x = or t = 0, S2M + Si/i = 0. 

Here the annihilator bundle over the front face is given with respect to both the 
coordinate system (|2.19l) and (|2.20[) , 

Thus, the two subbundles over the old boundary hypersurfaces meet transversally 
over the intersection, up to the corner, as claimed and so can be blown up in either 
order. In the complement of the lifts of the old boundaries under these two blow ups, 
the variables fJ>/s\ and /// s 2 become legitimate; in terms of these the subbundle 
over the front face becomes smooth up to, and with a product decomposition at, 
all its boundaries. Thus, it too can be blown up. That the result is a (painful) 
reconstruction of the b-cotangent bundle of the blown up manifold Mq follows 
directly from the construction. 

It remains to consider the inclusion into b,ff TM 2 . The only changes are at the 
front face, namely the third line of (|2.23p becomes 
(2.26) 

lsd s + v(sx'd y - sd Y ) + l'(x'd x > - sd s - Yd Y ) + v'dy = at x' = 0, s > 
<=*l = l', v' = sv + I'Y. 
Correspondingly the third line of (|2.25[) becomes 

x' = 0, A + A' + f/Y = 0, fi + sn' = 
( 2 - 27 ) or t = 0, s 2/ u + sim' = 0, s 2 (A + A') tfj Z i = °- 

The rest of the argument is unchanged, except that the conclusion is that h S T*MQ 
is being reconstructed. □ 

Note that for any manifold with corners, X, the b-cotangent bundle of any 
boundary hypersurface H (or indeed any boundary face) includes naturally as a 
subbundle h T*H ^ h T^X. 

Lemma 2.6. The Hamilton vector field of gs lifts from either the left or the right 
factor of M in M 2 to a smooth vector field, tangent to the boundary hypersurfaces, 
on h T*Ml, as well as on b ' T*Mq, still denoted by H^, resp. H^. Moreover, 
Hn, = PlV l , = prV r , where Vl, resp. Vr are smooth vector fields tangent to 
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all hypersurfaces except the respective cotangent bundles over the left, resp. right, 
boundaries, to which they are transversal, and where pl and pn are defining func- 
tions of the respective cotangent bundles over these boundaries. 

Proof. Inserting the explicit form of the Euclidean metric, the Hamilton vector field 
in (I2.15[) becomes 

(2.28) °H po = X(xd x + pd^) - \p\ 2 d x + xp ■ d y . 

Consider the lift of this vector field to the product, °T*M x °T*M, from left and 
right, and then under the blow up of the diagonal in the boundary. In the coordinate 
systems (j2~19|) and 

°H£ = X(sd s + lldfj) - \p\ 2 d X + sfi(x'd y - dy) 

°H« = \'{x'd x i - S8 S - Ydy + - \p!\ 2 d X , + p'8y 

°H£ = \( Sl d Sl + fidfj) - \p\ 2 8 x + siJ2 N{^ Vi - d Zj ) 

(2.29) +sipi(td Vl - td t + sid Sl + s 2 d S2 + ^ Z i d Zj) 
°H£ = \'{s 2 d S2 + - \p'\ 2 dy +S2^2 Hi 8 ** 

+s 2 p'i{td t - s 1 d Sl - s 2 d S2 - ^ z i d Zj)- 

J>2 

Note that the bundle itself is just pulled back here, so only the base variables are 
changed. 

Next we carry out the blow ups of Lemma l2.5l The centers of blow up are given 
explicitly, in local coordinates, in (|2.25[) . with the third line replaced by (|2.27[) in 
the case of l 2 ^ s . We are only interested in the behaviour of the lifts of the vector 
fields in (|2.29[) near the front faces introduced in the blow ups. 

Consider tg b first. For the first two cases there are two blow-ups, first of /z = 
in s = and then of p + sp' = in x' = 0. Thus, near the front face of the first 
blow up, the p variables are replaced by p = p/s and then the center of the second 
blow up is p + p' = 0, x' = 0. Thus, near the front face of the second blow up we 
can use as coordinates s, x' , p' and v = [p + p')/x' , i.e. substitute p = —p' + x'v. 
In the coordinate patch (|2.29[) the lifts under the first blow up are 

°H£ ^ \sd s - S 2 \p\ 2 3 X + S 2 p{x'd y - dy) 
U " >1 " °H« .— > \'{x'd x , - sd s - Ydy + pdft + p'd^) - \p'\ 2 d X > + p'dy. 

Thus under the second blow up, the left Hamilton vector field lifts to 

(2.31) °H£ o .— > sT, T = \d s - s\p\ 2 d x + sp(x'd y - dy) 

where T is transversal to the boundary s — where A ^ 0. 

A similar computation near the corner shows the lifts of the two Hamilton vector 
fields under blow up the fibrations of s\ = and s 2 = in terms of the new 
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coordinates ft, = /s/si and p! = p! jsi to be 

°H^ o =\ Sl d Sl - sl\il\ 2 d x + slJ2^(td w - d Z3 ) 

J>2 



(2.32) 



+s\jli(tdy 1 - td t + s 1 d Sl - fidf, + s 2 d S2 - ji'd^ + ^ z j d z j ), 

J>2 

»u« _\'„_« „2|r-,'|2c,. . , JVr,'. 



J H£ =A' S2 5 S2 - sll/i'I 2 A , + s 2 ^ V 3 dz 3 

i>2 

+s|/i' 1 (t9t - sid Sl + - s 2 d S2 + #d w - ^ Zjd Zj ). 

j>2 

The final blow up is that of i = 0, /2 + /2' = 0, near the front face of this blow-up 
replacing (t,p,,jl') by (t,ji,i>), v — {ji + fi')/t, as valid coordinates (leaving the 
others unaffected). Then the vector fields above become 

°H£ = Sl f \ °H« = aa f« 

f L =A0 Sl - Sl |/l| 2 a A + si^Mi(*^ - 3z 3 ) 

i>2 



(2.33) 



-sijli{tdy 1 - td t + sid Sl - //% + s 2 d S2 + ^ z i®Zj)i 



J>2 



f R =\'d S2 - s 2 \ii'\ 2 dy + s 2 J2^9 Zj 

+s 2 p' 1 (td t - si9 Sl + jldf, - s 2 d S2 - ^ z jdzj)- 

J">2 

Thus both left and right Hamilton vector fields are transversal to the respective 
boundaries after a vanishing factor is removed, provided A, A' 0. 

The final step is to show that the same arguments apply to the perturbed metric. 
First consider the lift, from left and right, of the perturbation to the Hamilton vector 
field arising from the perturbation of the metric. By assumption, the perturbation 
H is a 2-cotensor which is smooth up to the boundary. Thus, as a perturbation 
of the dual metric function on °T*M it vanishes quadratically at the boundary. 
In local coordinates near a boundary point it follows that the perturbation of the 
differential of the metric function is of the form 

(2.34) dp — dpo = x 2 (a h bdy + cdX + ed/j,) 

x 

From (|2.14[) it follows that the perturbation of the Hamilton vector field is of the 
form 

(2.35) Hp - H po = x 2 (a'xd x + b'xd y + c'dx + e'fy) 

on °T*M. Lifted from the right or left factors to the product and then under the 
blow-up of the diagonal to Mq it follows that in the coordinate systems (|2.19|) and 
(|2.20p . the perturbations are of the form 

3 Hp - H£ o = s 2 {x') 2 V L , H* - H« = (x'fV R , V L , V R e V b , 

( "" >(,! H£ - Hp Q = s\t 2 W L , H R - H R = s 2 t 2 W R , W R , W R G V b 

where Vb denotes the space of smooth vector fields tangent to all boundaries. Since 
the are lifted from the right and left factors, V L and V R are necessarily tangent 
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to the annihilator submanifolds of the right and left boundaries. It follows that 
the vector fields sx'V L and x'V R are tangent to both fibrations above a coordi- 
nate patch as in (gUgl and sxtV L and s 2 tV R are tangent to all three annihilator 
submanifolds above a coordinate patch (|2.20[) . Thus, after the blow ups which 
reconstruct b T*Mg, the perturbations lift to be of the form 

(2-37) H£ - H£ o = PLPs U L , H R - H R = p RPs U R 

where U L and U R are smooth vector fields on h T*M^. 

From (|2.37[) it follows that the transversality properties in (|2.31l) and (|2.33p 
persist. 

Now consider ff . First, (|2.30[) is unchanged, since the annihilators on the 'old' 
boundary faces are the same in this case. In particular, we still have p = fi/s as one 
of our coordinates after the first blow up; the center of the second blow up is then 
x' = 0, A + A' + // • Y = 0, p + // = 0. Thus, near the front face of the second blow 
up we can use as coordinates s, x', p! and cr = (A + A' + // • Y)/x', v = (p + //)/#', 
i.e. substitute p = —fjf + x'v, i.e. p = —p!s + x' sv, and A = —A' — // • Y + x'er. 
Thus under the second blow up, the left Hamilton vector field lifts to 

(2.38) °H£ o .— ► sT', T' = Xd s + p ■ {-vd a + x'd y - dy), 

so T' is transversal to the boundary s — where A 7^ 0. 

In the other coordinate chart, again, (|2.32[) is unchanged since the annihilators 
on the 'old' boundary faces are the same. The final blow up is that of 

t = 0, p. + A' = 0, A + A' + p[ + J2 frjZj = 0, 

i>2 

near the front face of this blow-up replacing (t, p, p! , A, A') by (t, p, u, A, a), 

v={p + p')/t, a = (A + A' + p[ + P>jZj)/t, 

i>2 

as valid coordinates (leaving the others unaffected). Then the vector fields above 
become 

f L =Xd Sl - Sl \p\ 2 d x + ax^piAtdy, - d Z] - i/jd a ) 

J>2 

+sipi(td Vl - td t + s\d si + ad„ - vdp + s 2 d s . 2 - v\d a + Zjd Zj ), 

(2.39) ~ j>2 

f R =X'd S2 +s 2 Y / ^dz 1 

J>2 

+s 2 p[(td t - s 1 d Sl + pdp, - s 2 d S2 - ada - 2J Zjd Zj )- 

j>2 

Again, both left and right Hamilton vector fields are transversal to the respective 
boundaries after a vanishing factor is removed, provided A, A' ^ 0. The rest of the 
argument proceeds as above. □ 

Proposition 2.7. The differential of the distance function for the perturbed metric 
gs, for sufficiently small 5, on M — B™ defines a global smooth Lagrangian subman- 
ifold Ag, of b T*(M Xo-M), which is a smooth section outside the lifted diagonal and 
which lies in b T* ff over the front face, ff(Mg) and in consequence there is a unique 
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geodesic between any two points of (B 3 )°, no conjugate points and (|2.3|) remains 
valid for dist^z, z'). 

Proof. For the unperturbed metric this already follows from Lemma 12.11 We first 
reprove this result by integrating the Hamilton vector fields and then examine the 
effect of the metric perturbation. Thus we first consider the lift of the Hamilton 
vector field of the hyperbolic distance function from °T*M, from either the left of 
the right, to b < ff T*M 2 , using the preceding lemma. 

Although the global regularity of the Lagrangian which is the graph of the dif- 
ferential of the distance is already known from the explicit formula in this case, 
note that it also follows from the form of these two vector fields. The initial mani- 
fold, the unit conormal bundle to the diagonal, becomes near the corner of M 2 the 
variety of those 

(2.40) £_(dx - dx') + rj{dy - dy') such that £ 2 + M 2 = — !-, x = x' > 0. y = y . 

2x z 

In the blown up manifold M5 2 the closure is smooth in b T*Mg and is the bundle 
over the lifted diagonal given in terms of local coordinates (|2.19[) by 

ds 1 

(2.41) A— + fxdY, A 2 + | M | 2 = -, s = 1, Y = 0; 

s 2 

the analogous statement also holds in h ' T*M{j where one has the 'same' expression. 
Using the Hamilton flow in b ' ff T*MQ, we deduce that the flow out is a global smooth 
submanifold, where smoothness includes up to all boundaries, of b ' ff T*Mg , and is 
also globally a graph away from the lifted diagonal, as follows from the explicit 
form of the vector fields. Note that over the interior of the front face, b,ff T*M ( 2 is 
just the standard cotangent bundle, so smoothness of the distance up to the front 
face follows. Over the left and right boundaries the Lagrangian lies in A = 1 and 
A' = 1 so the form (|2.3| of the distance follows. 

The analogous conclusion can also be obtained by using the flow in b T*Mg. In 
this setting, we need that over the front face (12.41[) is contained in the image of 
b T* ff to which both lifted vector fields are tangent. Thus it follows that the flow out 
is a global smooth submanifold, where smoothness includes up to all boundaries, of 
h T*M§ which is contained in the image of h T* ff over S(M§ ). Again, it is globally a 
graph away from the lifted diagonal, as follows from the explicit form of the vector 
fields. Over the left and right boundaries the Lagrangian lies in A = 1 and A' = 1 
so the form (|2.3p of the distance follows. 

For small S these perturbation of the Hamilton vector fields are also small in 
supremum norm and have the same tangency properties at the boundaries used 
above to rederive (|2.3p . from which the Proposition follows. □ 

3. The semiclassical double space 

In this section we construct the semiclassical double space, Mo.ft, which will be 
the locus of our parametrix construction. To motivate the construction, we recall 
that Mazzeo and the first author [TB] have analyzed the resolvent R(h, a), defined in 
equation (fTT|) . for a/h 6 C, though the construction is not uniform as \a/h\ — > oo. 
They achieved this by constructing the Schwartz kernel of the parametrix G(h, a) 
to R{h 1 a) as a conormal distribution defined on the manifold 

M =B" +1 x B" +1 



20 



RICHARD MELROSE, ANTONIO SA BARRETO AND ANDRAS VASY 



defined in ()2.2[) . see also figure Q] with meromorphic dependence on ajh. 

The manifold B™ +1 x B™ +1 is a In + 2 dimensional manifold with corners. It 
contains two boundary components of codimension one, denoted as in [16] by 

d[{M n+1 x B Il+1 ) = <9B" +1 x B™ +1 and <9[(B™ +1 x B" +1 ) = B™ +1 x ,9B™ +1 , 

which have a common boundary d 2 {M n+1 x B" +1 ) = dM n+1 x <9B n+1 . The lift of 
d[(B n+1 x B ,l+1 ) to M , which is the closure of 

(3- 1 (dl(M n+1 x B" +1 ) \<9 2 (B n+1 x B" +1 )), 

with (3 : Ma —> B n+1 x B™ +1 the blow-down map, will be called the right face and 
denoted by K. Similarly, the lift of <9j (B™ +1 x B" +1 ) will be called the left face and 
denoted by C. The lift of Diag ri(?2(B™ +1 x B™ +1 ), which is its inverse image under 
j3, will be called the front face J 7 , see figure [TJ 

We briefly recall the definitions of their classes of pseudodifferential operators, 
and refer the reader to [16 for full details. First they define the class ^^(B™" 1 " 1 ) 
which consists of those pseudodifferential operators of order m whose Schwartz 
kernels lift under the blow-down map j3 defined in (|2.2[) to a distribution which 
is conormal (of order m) to the lifted diagonal and vanish to infinite order at all 
faces, with the exception of the front face, up to which it is C°° (with values in 
conormal distributions). Here, and elsewhere in the paper, we trivialized the right 
density bundle using a zero-density; we conveniently fix this as |d<7j(,z')|. Thus, 
the Schwartz kernel of A £ *™(B n+1 ) is K A (z, z')\dg 5 {z% with K A as described 
above, so in particular is C°° up to the front face. 

It then becomes necessary to introduce another class of operators whose kernels 
are singular at the right and left faces. This class will be denoted by ty Q n ' a ' b (M n+1 ), 
a,beC. An operator P £ V n ' a ' b (M n+1 ) if it can be written as a sum P = P x + P 2 , 
where Pi £ *i>™(X) and the Schwartz kernel Kp 2 \dgs{z')\ of the operator P 2 is such 
that Kp 2 lifts under j3 to a conormal distribution which is smooth up to the front 
face, and which satisfies the following conormal regularity with respect to the other 
faces 

(3.1) VbP*Kp 2 £ p a L p b R L oc> (B n+1 x B" +1 ), V k e N, 

where Vfc denotes the space of vector fields on Mq which are tangent to the right 
and left faces. 

Next we define the semiclassical blow-up of 

B n+i x M n+i x j 0;1 ^ 

and the corresponding classes of pseudodifferential operators associated with it 
that will be used in the construction of the parametrix. The semiclassical double 
space is constructed in two steps. First, as in [TBJ, we blow-up the intersection 
of the diagonal Diag x[0, 1) with dM n+1 x dM n+1 x [0, 1). Then we blow-up the 
intersection of the lifted diagonal times [0, 1) with {h = 0}. We define the manifold 
with corners 

(3.2) M ,r= [B n+1 xB n+1 x [0,l] h ;dDiagx[0,l);Diag x{0}]. 
See Figures [5] and [3J We will denote the blow-down map 

(3.3) n : M ,h — ► B" +1 x B" +1 x [0, 1). 
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Figure 2. The stretched product B n+1 x B™ +1 x [0, 1). 




Figure 3. The semiclassical blown- up space Mq^ obtained by 
blowing-up B n+1 x B ,l+1 x [0, 1) along Diag n x '{h = 0}. 

As above, we can define the right and left semiclassical faces as the lift of 

8[ (B' i+1 x B" +1 x [0, 1)) = <9B" +1 x B n+1 x [0, 1) and 

d{ (B n+1 x B" +1 x [0, 1))) = B n+1 x dM n+1 x [0, 1), 

by the blow-down map [3%. We will denote the lift of the diagonal under the map 
/9fi by Diag,>, i.e. 

Diag„ = the closure of fi h 1 (Diag x(0, 1) \ (Diag n(<9B™ +1 x dW l+1 ))) . 

The lift of d Diag x [0, 1) will be called the zero front face T, while the boundary 
obtained by the blow-up of Diag x[0, 1] along Diag x{0} will be called the semi- 
classical front face S. The face which is obtained by the lift of B™ +1 x B™ +1 x {0} 
is the semiclassical face, and will be denoted by A. 

We wish to find a parametrix such that P(h, a) acting on the left produces the 
identity plus an error which vanishes to high enough order on the right and left 
faces, 1Z and £, and to infinite order at the zero-front face J 7 , the semiclassical front 
face S and the semiclassical face A. Thus the error term is bounded as an operator 
acting between weighted L 2 (B n+1 ) spaces and its norm goes to zero as h X 0. 

As in [TB] we define the class of semiclassical pseudodifferential operators in two 
steps. First we define the space P e ^Q l h (M> n+1 ) which consists of operators whose 
kernel 

K P {z,z',h)\dg s (z')\ 

lifts to a conormal distribution of order m to the lifted diagonal and vanishes to 
infinite order at all faces, except the zero front face, up to which it is C°° (with 
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values in conormal distributions) and the semiclassical front face, up to which it is 
fa-n-iQoc ( w ^h values in conormal distributions). We then define the space 

(3.4) 

A-''"'(.\/„,,; = {K e L°°(M 0tn ) : V?K e p^VflPs™ - ^ 00 ^), m g N}, 

where Vb denotes the Lie algebra of vector fields which are tangent to C, A and 
1Z. Again, as in [16], we define the space v[/™^ a,b ' c (B n+1 ) as the operators P which 
can be expressed in the form P = Pi + P^, with Pi g ^ Q n n (M n+1 ) and the kernel 
Kp 2 \dg s {z')\ of P 2 is such p* h Kp 2 g JC a ' b - c (M . h ). 

4. A SEMICLASSICAL PARAMETRIX FOR THE RESOLVENT IN DIMENSION THREE 

In this section we construct a parametrix for the resolvent R(h, a) defined in 
(|1T|) in dimension three. We do this case separately because it is much simpler 
than in the general case, and one does not have to perform the semiclassical blow- 
up. Besides, we will also need part of this construction in the general case. More 
precisely, if n + 1 = 3, and the metric g satisfies the hypotheses of Proposition 
12 . 71 we will use Hadamard's method to construct the leading asymptotic term of 
the parametrix G(tr, h) at the diagonal, and the top two terms of the semiclassical 
asymptotics. Our construction takes place on (B 3 Xo IB 3 ) x [0, l)h, instead of its 
semiclassical blow up, i.e. the blow up of the zero-diagonal at h — 0, as described 
above. This is made possible by a coincidence, namely that in three dimensions, 
apart from an explicit exponential factor, the leading term in the asymptotics lives 
on B 3 x B 3 x [0, 1). However, to obtain further terms in the asymptotics would 
require the semiclassical blow up, as it will be working in higher dimensions. For 
example, this method only give bounds for the resolvent of width 1, while the more 
general construction gives the bounds on any strip. 

We recall that in three dimensions the resolvent of the Laplacian in hyperbolic 
space, A go , Ro(a) = (A 9o — a 2 — 1) _1 has a holomorphic continuation to C as an 
operator from functions vanishing to infinite order at <9B 3 to distributions in B 3 . 
The Schwartz kernel of Rq (<t) is given by 

iar 

(4-1) Ro(a,z,z')= , 

47r sinn r$ 

where ro = Tq(z, z') is the geodesic distance between z and z' with respect to the 
metric go, see for example [16] , 

Since there are no conjugate points for the geodesic flow of g, for each z' g 
(B 3 ) , the exponential map for the metric g, exp g : IV(B 3 ) — > (B 3 ) , is a global 
diffeomorphism. Let (r, 9) be geodesic normal coordinates for g which are valid 
in (B 3 )° \ {z 1 }; r(z,z') =: d(z,z') is the distance function for the metric g. Since 
r(z, z') is globally defined, g is a small perturbation of go and the kernel of i?o( cr ) 
is given by (|4.1I) . it is reasonable to seek a parametrix of R(h, a) which has kernel 
of the form 

(4.2) G(h, a, z, z') = e- l ^ r h- 2 U(h, a, z, z% 

with U properly chosen. 

We now reinterpret this as a semiclassical Lagrangian distribution to relate it to 
the results of Section [21 Thus, —err = —crd(z,z') is the phase function for semi- 
classical distributions corresponding to the backward left flow-out of the conormal 
bundle of the diagonal inside the characteristic set of 2p c — a 2 . This flowout is the 
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Figure 4. The manifold B™ +1 Xi B" +1 . 

same as the forward right flow-out of the conormal bundle of the diagonal, and is 
also the dilated version, by a factor of a, in the fibers of the cotangent bundle, of 
the flow-out of in the characteristic set of 2p e — 1 , which we described in detail in 
Section H 

In view of the results of Section^ for the characteristic set of 2p e — 1 (the general 
case of 2p c — a 2 simply gives an overall additional factor of a due to the dilation), 
the lift p* d r of d r to B 3 x B 3 satisfies 

(4.3) p*d r = n e *v e L , 

and thus is a C°° vector field on (B 3 XoB 3 ) \Diag which is tangent to all boundary 
faces, and at the left face L = B, 

(4.4) d r = -R L + W L , 

where is the radial vector field corresponding to the left face, and Wl S 
p L V b ((B 3 x B 3 )\Diag ). 

It is convenient to blow up Diag and lift (3*d r further to this space. We thus 
define B 3 xi B 3 to be the manifold obtained from B 3 Xo B 3 by blowing-up the 
diagonal Diag as shown in Fig. |4] Let (3d : B 3 XiB 3 i — > B 3 x B 3 denote the blow- 
down map and let fto (3 = /3od- The vector field P^dr is transversal to the new face 
introduced by blowing up the diagonal; it still satisfies the lifted analogue of (14. 4[) . 
Moreover, integral curves of 0Q d d r hit the left face L away from its intersection with 
the right face R in finite time. 

In coordinates (r, 6) the metric g is given by 

(4.5) g = dr 2 + H(r,9,d6), 

where H (r, 9, d9) is a C°° 1-parameter family of metrics on § 2 . The Laplacian with 
respect to g in these coordinates is given by 

A g = -d 2 r -Vdr + A H , V = -^drdgfi), 

\g\ 2 

where \g\ i is the volume element of the metric g and A# is the Laplacian with 
respect to H on S 2 . The volume element \g\^ has the following expansion as r ! 0, 

(4.6) \g\^{r,9) = r 2 {l + r 2 gi {r,9)), 
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see for example page 144 of j7]- So 

(4.7) A g = -d 2 r - + rA)d r + A H 



We want U in (14. 2|) to be of the form 

(4.8) U(h,a,z,z') = U (a,z,z') + hUi(a,z,z'), 
and so 

(h 2 (A g + x 2 W - 1) - a 2 ) e-^ r h- 2 U 

(4.9) = ^ r ((A ff + * 2 VF - 1)U + 2^| 5 |-^ r (| 5 |*?7o) 

+ 2i<r| 5 |-26l r .(|g|iC/i) + fc(A fl + x 2 T^ - l)C7i) . 
Here the leading term in h as h — >■ both overall, and as far as [To is concerned, 



is 



(4.10) 2i£| ff |-*0 r (| ff |*tfo), 
and the leading term as far as U± is concerned is 
(4.H) 2ia\g\-id r (\g\iu 1 ). 

In the interpretation as semiclassical Lagrangian distributions, these are both the 
differential operators arising in the transport equations. For Hamiltonians given 
by Ricmannian metrics, these operators in the interior are well-known to be Lie 
derivatives with respect to the Hamilton vector field, when interpreted as acting on 
half-densities. This can also be read off directly from (|4. 10|) - (|4. 1 1|) . with \g\* being 
the half-density conversion factor, a is due to working at energy a (rather than 1), 
and the factor of 2 is due to the symbol of the Laplacian being 2p e . 

To get rid of the term in h~ l we solve the 0th transport equation, i.e. we impose 

d r (\g\iu ) = 0, 
and we choose Uo(r,0) — ■r-\g(r, From (|4.6|) we have 

(4.12) \g\-i(r,d) = r _1 (l + r 2 g 2 {r, 9)) near r = 0. 
Therefore, near r — 0, 

(4.13) A ff i-| 9 (r,0)r* = S(z, z<) + ^-Ar' 1 + ^A g (rg 2 ). 

This only occurs in three dimensions, and makes this construction easier than in 
the general case. In higher dimensions the power or r in f|4. 1 2[) is r~ ~s , and does not 
coincide with the power of r of the fundamental solution of the Laplacian, which, 
in dimension n + 1, is r 1- ", so one does not get the delta function in (|4.13l) . 

To get rid of the term independent of h in (|4.9j) in r > we solve the first 
transport equation, 

2ia\g\-id r (\g\iu 1 ) + (A g + x 2 W - 1)U = in r > 0, 
Ut = at r = 0. 

So 

(4.14) U 1 (r,0) = --^—\g(r,9)\-^[ \g\i(s, 9) (A g + x 2 W - l) \g\-*(s, 6) ds. 
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Since |<?| 4 is C°° up to r = 0, and vanishes at r = 0, it follows from (14. 13)) that 
| <7 1 4 A g |g|~4 is C°° up to r = 0. In particular the integrand in (|4.14j) is smooth up 
to r = 0. With these choices of Uo and Ui we obtain 

(h 2 (A g + x 2 W - 1) - a 2 ) e-^ r h- 2 U(h,a,z,z') 

' 4 ' ' ' ' = 5(z, z') + he-^ r (A g + x 2 W - l)U x {(j, z, z') 

This gives, in principle, a parametrix G(h,a, z, z') = e~ l % r h~ 2 U(h,a, z, z') in 
the interior of B 3 x B 3 in the two senses that the diagonal singularity of R(h, a) is 
solved away to leading order, which in view of the ellipticity of the operator means 
that the error 

E(h, a) = (h 2 {A g + x 2 W - 1) - a 2 ) G(h, a) - Id 

is a semiclassical pseudodifferential operator of order —1 (in a large calculus, i.e. 
with non-infinite order vanishing off the semiclassical front face, which did not even 
appear in our calculations), and the top two terms of the semiclassical parametrix 
of (h 2 {A g + x 2 W - 1) - cr 2 ) -1 are also found. 

In fact, our parametrix is better than this. To understand the behavior of G and 
the remainder 

(4.16) E(h, a, z, z') = he-^ r (A g + x 2 W - l)U x 

near the boundary of B 3 x B 3 , we need to analyze the behavior of Uq and U\ 
at the left and right boundary faces. We will do the computations for arbitrary 
dimensions, since we will need some of these estimates in the general case, but in 
this special situation we have n = 2. 

We start by noting that the asymptotics of Uo and U\ follow from the transport 
equation which they satisfy. Indeed, much like we analyzed the flow-out of the 
conormal bundle of the diagonal, we show now that 

(4.17) U e p^pl^p^C 00 ^ 1 X!B ,l+1 ), Ui G R 2 pl /2 pl /2 C°°(M n+1 xjl^ 1 ); 

here po is the defining function of the front face of the blow-up creating B™ +1 Xj 
B™ +1 . Note that we have already shown this claim near this front face; the main 
content of the statement is the precise behavior as pl,Pr 0. 

We start with Uq. First, the conclusion away from the right face, pn = 0, follows 
immediately from (|4.4I) since integral curves emanating from the lifted diagonal hit 
this region at finite time, and solutions of the Lie derivative equation have this form 
near the boundary. To have the analogous conclusion away from the left face, we 
remark that solutions of the left transport equation automatically solve the right 
transport equation; one can then argue by symmetry, or note directly that as — d r > is 
the radial vector field at the right face, modulo an element of /9^Vb(B n+1 XoB™ +1 ), 
and as integral curves of d r > emanating from the lifted diagonal hit this region at 
finite time, and solutions of the Lie derivative equation have this form near the 
boundary. It remains to treat the corner where both pl = and pr = 0. The 
conclusion here now follows immediately as integral curves of d r reach this corner 
in finite time from a punctured neighborhood of this corner, and in this punctured 
neighborhood we already have the desired regularity. This proves (|4.17[) for Uq. 

To treat Ui, it suffices to prove that 

(4.18) (A 9 + x 2 W - 1)U G i? 2 p" /2+2 p™ /2 C 00 (B™ +1 x B" +1 \ Diag ), 
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for then 

(4.19) Ui G i? 2 p" /2 p™ /2 C 00 (B"+ 1 x B" +1 \ Diag ), 

by the same arguments as those giving the asymptotics of Uo, but now applied to 
the inhomogeneous transport equation. 
On the other hand, (|4. 18|) follows from 

/3*tt2(A s + x 2 W - 1) G Diffg(B n+1 x B" +1 ) 

with 

f3*TT* L ((A g + x 2 W-l)-(A go -1)) ei? 2 p|Diff 2 (B" +1 x B" +1 ); 

here ttt, is added to emphasize the lift is that of the differential operator act- 
ing on the left factor; lifting the operator on the right factor results in an 'error' 
i? 2 p|.Diff 2 (B™ +1 x B" +1 ). These two in turn follow immediately from the form of 
the metric, namely g - g e G x 2 C°°{X; °T*X <g> °T*X). 

This completes the proof of (|4. 1 T[> . and also yields that, with n + 1 = 3, 

(4.20) [3* {{A g + x 2 W - l)?7i) G i? 2 ^ /2+2 p^ /2 C°°(B n + 1 x B n+1 \ Diag ). 

Therefore, in the case n + 1 — 3, we have proved the following 

Theorem 4.1. There exists a pseudodifferential operator, G(h,a), a ^ 0, whose 
kernel is of the form 

G(h, a, z, z') = e~ l Tth~ 2 (U (h, a, z, z') + hUi(h, a, z, z')) 

with Uq and U± satisfying (14. 17)) and such that the error E(h, a) = P(h, a)G(h, a) — 
Id is given by (|4.16[) and satisfies (|4.20|) . 

5. The structure of the semiclassical resolvent 

In this section we construct the general right semiclassical parametrix G(h, a) 
for the resolvent. We will prove the following 

Theorem 5.1. There exists a pseudodifferential operator G(h,a), a =fi 0, such that 
its kernel is of the form 

(5.1) G(h,a,z,z') = e- l ^U(h,<T,z,z'), U E (M n+1 ). 

and that, using the notation of section 

P{h,a)G(h,a)-ld G P^Ps%. ;~'°° ,00 '* +i * (B" +1 ). 

Now, r 2 is a C°° function on the zero double space away from the left and 
right faces and in a quadratic sense it defines the diagonal non-degenerately; r has 
an additional singularity at the zero-diagonal. Correspondingly (^-j is C°° on 
Mo,n away from C, 1Z and A and defines the lifted diagonal non-degenerately in a 
quadratic sense; has an additional singularity at the zero diagonal. In particular, 
e~ l Tr is C°° on Mo.ft away from C, 1Z, A and the lifted diagonal, Diag R ; at Diag R 
it has the form of 1 plus a continuous conormal function vanishing there. Thus, 
its presence in any compact subset of Mq,s \(£U1ZU A) is not only artificial, but 
introduces an irrelevant singularity at Diag^, so it is better to think of G as 

G(h, a, z, z') = G' + e- iS ?U'(h, a, z, z'), 
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where G' is supported near Diag s (i.e. its support intersects only the boundary 
hypersurfaces S and but not the other boundary hypersurfaces) , and U' vanishes 
near Diag s . 

Indeed, the first step in the construction of G is to construct a piece of G', 
namely to find Go G such that Eq = P(h, a)Go(h, a) — Id has no singularity 
at the lifted diagonal, Diag^, with Go (hence Eq) supported in a neighborhood of 
Diag ft in Mo,b that only intersects the boundary of Mq^ at iS and J- . Thus, 

(5.2) P(h, a)G - Id = E , with E G %™(B n+1 ). 

Since P(h,a) is elliptic in the interior of B" +1 , the construction of Go in the 
interior follows from the standard Hadamard parametrix construction. We want to 
do this construction uniformly up to the zero front face and the semiclassical front 
face in the blown-up manifold Mofi. We notice that the lifted diagonal intersects 
the boundary of Mo t % transversally at the zero and semiclassical front faces. Since 
Diag R intersects the boundary of Mq^ transversally at S and F, we proceed as 
in [15], and extend Diag R across the boundary of Mq^, and want to extend the 
Hadamard parametrix construction across the boundary as well. To do that we 
have to make sure the operator P(h, a) lifts to be uniformly transversally elliptic 
at the lift of Diag^ up to the boundary of Mo,fc, i- e - it is elliptic on the conormal 
bundle of this lift, and thus it can be extended as a transversally (to the extension 
of the lifted diagonal) elliptic operator across the boundaries of Mq^. Note that up 
to S \ J 7 , this is the standard semiclassical elliptic parametrix construction, while 
up to T \ S, this is the first step in the conformally compact elliptic parametrix 
construction of Mazzeo and the first author [TB], so the claim here is that these 
constructions are compatible with each other and extend smoothly to the corner 
S H T near Diag^. 

To see the claimed ellipticity, and facilitate further calculations, we remark that 
one can choose a defining function of the boundary x such that the metric g can 
be written in the form 

dx 2 H e (x,ui) 

•9 = — H 2 ' 

x z x z 

where H e is one-parameter family of C°° metrics on S n . In these coordinates the 

operator P(h, a) is given by 

(5.3) 



P(h,a) = h 2 ~{xd x ) 2 + nxdn + x 2 A(x,uj)d x + x 2 A He ( x ^) + x 2 W - 




We then conjugate P(h, a) by x 2 , and we obtain 

(5.4) Q(h,a) = x~%P(h,a)x* = h 2 {-{xd x f + x 2 Ad x + x 2 A Hc +x 2 B) - a 2 , 

where B = — '^A + W. To analyze the lift of Q(h, a) under we work in projective 
coordinates for the blow-down map. We denote the coordinates on the left factor of 
jjn+i while the coordinates on the right factor will be denoted by (x 1 , lu'). 

Then we define projective coordinates 

(5.5) x' = P , x = Z, r = ^, 

x' x' 

which hold away from the left face. The front face is given by7 r ={/? = 0} and 
the lift of the diagonal is Diag = {X = l,Y = 0}. The lift of Q(h,a) under the 
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zero blow-down map (3 is equal to 

Qo{h,a)=j3*Q(h,a) = 
h 2 {-(Xd x ) 2 +X 2 Apd x +X 2 A HApX ^ +pY) (D Y ) -X 2 p 2 B) -a 2 . 

In this notation, the coefficients of A H / pXu / +pY ^ (Dy) depend on p, oj' and Y, but 
the derivatives are in Y. This operator is transversally elliptic in a neighborhood of 
{X = 1, Y = 0}, away from h = 0. 

The restriction of the lift of Qo(h, a) to the front face J- = {p = 0}, is given by 

NAQo(h,a)) = h 2 (-(Xd x ) 2 + X 2 A HAQ ^ } (D Y )) - a 2 . 

As in [IB] , Njr(Q (h, a)) will be called the normal operator of Qo(h, a) at the zero 
front face T . Notice that it can be identified with Laplacian with respect to the 
hyperbolic metric on the half-plane {X > 0,Y G R™} with metric X~ 2 (dX 2 + 
H e (0,oj')) conjugated by I>. 

Now we blow-up the intersection of Diag x [0, h) with h — 0. We define projective 
coordinates 

(5.6) ^ = ^> Y * = T 

The lift of Q(h,a) under the semiclassical blow-down map is given in these 
coordinates by 

Q h = /3* h Q(h,a) 

= -((1 + hX h )d Xh ) 2 + (1 + hX h ) 2 A He (D Yh ) 

- (1 + hX h ) 2 Aphd Xh + h 2 p 2 (l + hX h ) 2 B - a 2 , 

where H e = H t (p(\ + hXri),uj' + phYf,). This operator is transversally elliptic to 
{X h = 0, Y h=0 } near h = 0. 

The restriction of the lift of Qn to the semiclassical face, S = {h = 0} will be 
called the normal operator of Qn at the semiclassical face, and is equal to 

(5-7) N s (Qh) = ~d 2 Xn + A Ht(piUl) (D YR ) - a 2 . 

This is a family of differential operators on R^-^y^ depending on p and ui 1 ', and for 
each ui' and p fixed, Ns(Qfi) + a 2 is the Laplacian with respect to the metric 

S R = dX 2 +Y,He,ij{p,u)dY H! idY Htj , 

which is isometric to the Euclidean metric under a linear change of variables 
(XfijYfi) for fixed {p,uj'), and the change of variables can be done smoothly in 
(p, w'). Note that the fibers of the semiclassical blow-down map fin on S are given 
exactly by (p, uj') fixed. 

Therefore, the operator Q(h, <r) lifts under to an operator which is elliptic 
in a neighborhood of the lifted diagonal uniformly up to the zero front face and 
the semiclassical front face. Since the diagonal meets the two faces transversally, 
one can extend it to a neighborhood of T and S in the double of the manifold 
M§fi across S and and one can also extend the operator to be elliptic in 
that neighborhood. Now, using standard elliptic theory (or, put somewhat dif- 
ferently, the standard theory of conormal distributions to an embedded subman- 
ifold without boundary, in this case the extension of the diagonal), one can find 
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a Go G ^ \(M n+1 ) whose Schwartz kernel lifts to a distribution supported in a 
neighborhood of Diag^ such that 

(5.8) Q(h, <t)G -ld = E e %^(M n+1 ). 

Next we will remove the error at the semiclassical front face. We will find an 
operator 

d = G[ + e-^ r U[, G[ G %^(M n+1 ), U[ G ¥~~' 00, ~3~ 1,0 °(B n+1 ) 
such that Gi is supported near Diag^ while U[ is supported away from it, and 

Q(h,a)G 1 -E = E u 

(5.9) E x =E' l +e- i ^F 1 , E[ G pf* ,ft, *i G -".oo.-f-i.oo 

and with (3^ 1 Ke 1 supported away from £, 1Z, 

and Kej, resp. supported near, resp. away from, Diag^. In other words, the 
the error term E\ is such that the kernel of E[ vanishes to infinite order at all 
boundary faces (hence from now on we can regard it as trivial and ignore it), while 
the kernel of F\ lifts to a C°° function which is supported near S (and in particular 
vanishes to infinite order at the right and left faces), vanishes to infinite order at 
the semiclassical front face, and also vanishes to order § at the boundary face A. 

We will use the facts discussed above about the normal operator at the semi- 
classical face, Ns(Qn)- Notice that S, the semiclassical front face, is itself a C°° 
manifold with boundary which intersects the zero front face, J-, transversally, and 
therefore it can be extended across T. Similarly, the operator Ns{Qn) can be ex- 
tended to an elliptic operator across JF . We deduce from (|5.7|) that for each p and 
u>' fixed, and for Imcr < 0, the inverse of Ns(Qn) is essentially the resolvent of 
the Euclidean Laplacian at energy a 2 , pulled back by the linear change of variables 
corresponding to H e (p,u)'); for Imcr > we use the analytic continuation of the 
resolvent from Im a < 0. Here is where we need to make a choice corresponding 
to the analytic continuation of the resolvent of P(h,a) we wish to construct, i.e. 
whether we proceed from Imcr > or Imcr < 0; we need to make the corresponding 
choice for the Euclidean resolvent. 

Let i?o denote the analytic continuation of the inverse L 2 —> H 2 of the family 
(depending on p, u/) Ng^Qn) from Imcr < 0; it is thus (essentially, up to a linear 
change of coordinates, depending smoothly on p, uj') the analytic continuation of 
the resolvent of the Euclidean Laplacian. Since we are working with the analytic 
continuation of the resolvent, it is not automatic that one can solve away exponen- 
tially growing errors which arise in the construction below (i.e. that one can apply 
Rq iteratively to errors that arise), and thus it is convenient to make the follow- 
ing construction quite explicit order in h we are merely in the 'limiting absorption 
principle' regime (i.e. with real spectral parameter), thus the construction below is 
actually stronger than what is needed below. Moreover, from this point of view the 
construction can be interpreted as an extension of the semiclassical version of the 
intersecting Lagrangian construction of |18j extended to the 0-double space; from 
this perspective the method we present is very 'down to earth'. 

Via the use of a partition of unity, we may assume that there is a coordinate 
patch U in B" +1 (on which the coordinates are denoted by z) such that Eq is 
supported in ^^{U x U x [0, 1)). Note that coordinate charts of this form cover 
a neighborhood of S, so in particular E is in C°°(B n+1 x B n+1 x [0, 1)) outside 
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these charts, hence can already be regarded as part of the final error term and we 
can ignore these parts henceforth. 

Now, near S, /3^ 1 (t^ x U x [0, 1)) has a product structure °TU x [0, 5 ) = U x 
jjn+i x [0,<5 ), where °TU denotes the fiber-compactificd zero tangent bundle, and 
[0, Sq) corresponds to the boundary defining function ps- Indeed, the normal bundle 
of Diag in B" +1 x B' i+1 can be identified with °TB n+1 , via lifting 0-vector fields 
from on B™ +1 x B™ +1 via the left projection, which are transversal to Diag , 
hence the interior of the inward pointing spherical normal bundle of Diag x {0} in 
M = B ,l+1 x B ,l+1 x [0, 1) can be identified with °TB' i+1 , while the inward pointing 
spherical normal bundle itself can be identified with the radial compactification 
of °TB" +1 . However, it is fruitful to choose the identification in a particularly 
convenient form locally. Namely, away from <9Diag x{0}, coordinates z on U give 
coordinates 

z',Z% = (z - z')/h,p s 
near the interior of the front face S (here Z^ is the coordinate on the fiber of Tj/B™ +1 
over z'), while near <9Diag x{0}, (x 1 , u/, Xft, Y%, p$) (see (|5 . 6|) ) are coordinates 
near the interior of S (now {X%, Y%) are the coordinates on the fiber of °TjjM n+1 
over {x 1 ,oj')). To obtain coordinates valid near the corner, one simply needs to 
radially compactify the fibers of °Tj/B™ +1 , i.e. replace the linear coordinates Z^, 
resp. (Xfi, Yfr) by radially compactified versions such as \Zn\ and Z% = Z%/\Zfi\ E S n 
in the former case. 

Moreover, if a function, such as Eq, is supported away from A, then its support 
is compact in the interior of the fibers B n+1 of the fiber-compactified tangent space. 
Now, the interior of B' i+1 is a vector space, T p U, p G U, and in particular one can 
talk about fiberwise polynomials. Over compact subsets of the fibers, the boundary 
defining function ps is equivalent to h, and indeed we may choose boundary defining 
functions pa and ps such that 

h = paps- 

Note that pa is thus a boundary defining function of the compactified fibers of 
the tangent bundle; it is convenient to make a canonical choice using the metric 
<? e , which is an inner product on °T p U, hence a translation invariant metric on 
the fibers of TU , namely to make the defining function pA the reciprocal of the 
distance function from the zero section (i.e. the diagonal under the identification), 

smoothed out at the zero section. In particular, if U is a coordinate chart near 

_ i 

dM n+1 then pa = ((^i) 2 + l^ilf? ) 2 ■ This is indeed consistent with our previous 
calculations since 

PA = h((X-lf + \Y\ 2 y l2 
and therefore it is, away from Diag^, a defining function of the semiclassical face 
A. 

If v G C°° (Afo.fi), then expanding v in Taylor series around S up to order N, we 
have 

v = £ P k s v'k + v', v' k G C°°(T(7), v' G p N s +1 C°°{M a , h ). 

k<N 

In terms of the local coordinates valid near the corner S Pi A over an interior 
coordinate chart U, 

i /i z ~ z h i 
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v k is a C°° function of \z — z'\, \ z _ z z i\ , z ' ■ It is convenient to rewrite this as 
(5.10) 

v = h k \z-z'\- k v' k +v' = h k v k +v', v k e p A k C™(TU), v' g ^ +1 C°°(Mo, S ), 



k<N k<N 



for the reason that the vector fields D Zj are tangent to the fibers given by constant 
h, i.e. commute with multiplication by h. One can rewrite v completely analogously, 

v = Y hky K + v '> v " E P A k C°°(TU), v' G p N s +1 C°°{M ofi ), 

k<N 

for coordinate charts U at <9B" +1 . 

In addition, if a € C°°(B n+1 x B" +1 x [0, 1)), then expanding a in Taylor series 
around Diag x{0}, shows that for any N, modulo p$ +1 C°°(U x B™ +1 x [0, S )), it 
is of the form 

Y a atk {z'){z-z') a h k = Y, a a , k (z')Z^h k+ ^ 

\a\+k<N \a\+k<N 

= Y a a , k (z')Z%p k ^p k A , 

\a\+k<N 

where Z% = Z%/\Z n \\ a \ (except near the zero section) is C°° on B n+1 . While the 
last expression is the most geometric way of encoding the asymptotics at dM n+1 , it 
is helpful to take advantage of the stronger statement on the previous line, which 
shows that the coefficients are polynomials in the fibers, of degree < N. 

The vector fields hD z , resp. hDx and hDy, acting on a modified Taylor series 
as in (|5.10p . become Dz h , resp. Dx h and Dy h , i.e. act on the coefficients v k only 
(and on v' , of course) so we obtain that, modulo coefficients in p$ +1 C°°(U x B' i+1 x 
[0, Sq)), P(h,cr) lifts to a differential operator with polynomial coefficients on the 
fibers, depending smoothly on the base variables, i.e. an operator of the form 

Y a a , k ,^z')Z^h k+ ^D Zh , 

\a\+k<N,\fi\<2 

with an analogous expression in the (X^, Y^) variables. Here the leading term in h, 
corresponding to h° , is Ns{Qn) = &-gJx',uj') ~ ^ i i- e - we have 

P{h. / a)h m v m 

= h m (A gcix ,^ ) -<T 2 )v m + h m Y a a ^(z')Z%h k +^D Zr v m . 

0<\a\+k<N,\/3\<2 

Thus, one can iteratively solve away E$ as follows. Write 

E = h n_1 Y, h m Eo, m 
as in (|5.10[) . and note that each E _ m is compactly supported. Let 

so 

P(^, ( r)/ l m -"- 1 G 1:0 , m 
= h m E , m +h m - n - 1 Y a a . k . fi (z')Z^h k+ ^- n - 1 D Zh RoE , m , 

0<\a\ + k<N,\l3\<2 



32 



RICHARD MELROSE, ANTONIO SA BARRETO AND ANDRAS VASY 



and thus wc have replaced the error h rn Eo, m by an error of the form 

0<|a|+fe<JV,|/3|<2 
1<£<N \a\<£\P\<2 

= £ h m + e -^L e R E , m 
i<e<N 

which has the feature that not only does it vanish to (at least) one order higher in 
h, but the h m+l - ( ~ n+v > term is given by a differential operator Lg with polynomial 
coefficients of degree < I applied to R E m , with E m compactly supported. As 
the next lemma states, one can apply Rq to an expression of the form LgRoEo^m-, 
and thus iterate the construction. 

Lemma 5.2. Suppose Mj, j — 1,2, ... , N , are differential operators with polyno- 
mial coefficients of degree rrij on R^, +1 . Then 

R M 1 RoM 2 . . . RqM n R : C C °°(R" +1 ) -> C-°°(IR" +1 ) 

has an analytic extension from Imcr < to Recr > 1, Im cr £ K, and 

R M 1 R M 2 . ..RqMnRq ■■ C C °°(IR" +1 ) -> e- ia{w) (w)- n/2+N+m C°°(M n+1 ), 

with m = ^2 m j ■ 

Proof. For Imcr < 0, D Wk commutes with A — a 2 , hence with Rq, while commuting 
D Wk through a polynomial gives rise to a polynomial of lower order, so we can 
move all derivatives to the right, and also assume that Mj = w a ' J) , G N™ +1 , 
\aW)\ < mj, so we are reduced to examining the operator 

Rqw Row ...Row Ro- 
ll is convenient to work in the Fourier transform representation. Denoting the 
dual variable of w by £; for Imcr < 0, Ro is multiplication by (|£| 2 — o -2 ) -1 ! while 
w a ' J) is the operator (—D^)" 3 . Rewriting ^(Row " 1 ^ Row al2) . . . i? u> Q<N) Rof), the 
product rule thus gives an expression of the form 

E (ICI 2 - ^ 2 )- (Ar+1+|Q<1)|+ "' +|Q(N,|) Qa,^(C)(-^)^/ 

| / 3|<a< 1 )+. ..+«(") 

£ del 2 - ( T 2 )- (Ar+1+|Q<1,|+ - +|a(N)|) ^Q a ,A r ,,(^)^/, 

\/3\<aW + . ..+«<«) 

where Q a ,N,p is a polynomial in £. Since we are considering compactly supported /, 
the differential operator Q a .N.p{D w )w^ is harmless, and we only need to consider 
Ro +1+m applied to compactly supported functions. This can be further rewritten 
as a constant multiple of d^ +m Ro, so the well-known results for the analytic con- 
tinuation of the Euclidean resolvent yield the stated analytic continuation and the 
form of the result; see Proposition 1.1 of [17]. □ 

Applying the lemma iteratively, we construct 

G 1>m = e-™^ {Z h )- n / 2+N+m G' hm , G' hm G C°°(U x B" +1 ), 
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such that 



P(h,a) V h m ~^G hm ~ V h m -^E 0>m , 



m— m— 

where the series are understood as formal series (i.e. this is a statement of the 
equality of coefficients). Borel summing (z) m h m G[ m = Pg"G'x m , and obtaining 
G% G C°°(U x B n+1 x [0,5o)) as the result, which we may arrange to be supported 
where ps is small, we deduce that (|5.9|) holds. 

The next step is to remove the error at the semiclassical face A. We want to 
construct 

G 2 = e~^ r U 2 , U 2 ep^*oJ' f ^ f ~ 1,f , 

such that 



(5.11) 



P{h,a)G 2 - e-^Fx = E 2 , E 2 = e-^ r F 2l 
°a „ n / 2 „ n / 2 r oo rn n + 1 v„k™+ 1 v in -m — r,°°\ri~°°>% 



F 2 e h°°p% p n j{ Z C°°(B™ +1 x B™ +1 x [0, 1)) = pf^ 



h 



In other words, we want the error to vanish to infinite order at the semiclassical face 
A and at the semiclassical front face £, and to vanish to order ^ at the left and right 
faces; the infinite order vanishing at S means that S can be blown-down (i.e. does 
not need to be blown up), which, together with h being the joint defining function 
of A and 5, explains the equality of the indicated two spaces. This construction 
is almost identical to the one carried out in section |4] We begin by observing 
that the semiclassical face A consists of the stretched product B n+1 Xo B™ +1 with 
Diag blown-up, which is exactly the manifold B I1+1 Xi B n+1 defined in section 2) 
Moreover, as Fi vanishes to infinite order at S, see (|5.9I) . the latter can be blown 
down, i.e. Fx can be regarded as being of the form 

(5.12) 

Fx € /r n/2-1 C°°(B fl+1 x B" +1 x [0, 1)), with K Fl supported away from £, K, 
and vanishing to infinite order at Diag x [0, 1). 

Now, Fx has an asymptotic expansion at the boundary face h — of the form 

oo 

Fx ~ h-%- 1 h j Fx, j: Fx tj e C°°(B" +1 x B" +1 ), 

3=0 

Fx j vanishing to infinite order at Diag x [0, 1), supported near Diag x [0, 1). 



So we think of Fx as an element of B n+1 Xi B" +1 x [0, 1), where the blow-up 
■gm+i Xi fljn+i^ wag defied m section 21 see figure HJ with an expansion 

oo 

Fx-h-^J^VFiJ, 

3=0 

Fx j G C°°(B n+1 Xi B" +1 ), vanishing to infinite order at D. 

So we seek U 2 ~ h% ■ h?U 2 j with U 2 j vanishing to infinite order at D, such that 

P(h, a)e- l ^ r U 2 - e~^ r Fx = e^R, 
R e h°°C co (M n+1 xx B" +1 ), vanishing to infinite order at D. 
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Matching the coefficients of the expansions we get the following set of transport 
equations 

2ia\g\-*d r {\g\iu 2 , ) = -F h0 , if r > 0, 
U 2fi = at r = 0, 

and for j > 1, 

2ia\g\-id r (\g\iu 2tj ) = (A + x 2 W - y fe-i ~ F ltj , if r > 0, 
f7 2 .j = at r = 0. 

Notice that Fj_ is compactly supported and, as seen in equations (|4.17[) . U 2 ,o S 
pf p|c oc (B"+ 1 XiB"+ 1 ). Moreover, as in (jiTgj) . one gets that (A+a: 2 ^- ^)[/ 2 , G 
i? 2 pj|p£ +2 C°°(B n+1 Xi B' i+1 ), thus one can solve the transport equation for U 2 ,x, 
and gets that U 2 ,i € pf p%C°° {W l+1 x 1 B" +1 ). 

One obtains by using induction that U 2 ,j G pf p%C°° {W l+1 x 1 M n+1 ), and (A + 
x 2 W- ^-)U 2 ,j G i? 2 pf +2 p|c°°(B n+1 X!B ,l+1 ), for all j. Then one sums the series 
asymptotically using Borel's lemma. This gives U 2 and proves (|5.1ip . 

The last step in the parametrix construction is to remove the error at the zero 
front face. So far we have 

P(h,a) [G -G 1 + G 2 )-U = E 2 , 

with 

E 2 = e-^ r F 2 , F 2 G ^*o J' t+2 ^' f . 
Now we want to construct G3 such that 

(5.13) P(h,a)G 3 -E 2 = E 3 G e~** r pf pf%™ ,2+f *°°' ? (M n+1 ). 
We recall from Proposition 12. 71 that r = — log(pflpi) + F, F > 0. So, 

exp(-i^r) = (p i? p L ) i ^'exp(-z^F). 
h h 

Therefore the error term E 2 in (15 . 1 1[) satisfies 

E 2 = /3* h E 2 G pf^pr% 2 L +f +lf C°°(M ,,). 

We write 

00 

e 2 ~J2& E ^> E ^ e ^ PApl +ii p?* +i * c°°(n 

0=0 

and we want to construct G3 such that 

00 

/8ft(G 3 )~X;^G3j, 
i=o 

and that 

iV>(Qft)G 3 j = ^, G 3J G p| +lf C°°(F). 

The asymptotic behavior in p^ and pi follows from an application of Proposition 
6.15 of [16], and the fact that N^{Qn) can be identified with the Laplacian on 
the hyperbolic space. Now we just have to make sure, this does not destroy the 
asymptotics at the faces S and A. But one can follow exactly the same construction 
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we have used above, now restricted to the zero front face instead of Mo,n to construct 
G3J vanishing to infinite order at the faces S and A. 

This gives a parametrix, G = G - Gi + G 2 - G 3 G that 
satisfies 

(5.14) P(h, a)G-Id = R = E 3 + E' 3 e p'gf%*Q™' 2+ * +i *' 0o '* +i * (B™ +1 ). 
Since 

^3 6 ^°°*oT' 

and S3 is supported away from Diag s . The last step in the construction is to 
remove the error term at the left face and it will be done using the indicial operator 
as in section 7 of |16j . Since in the region near the left face is away from the 
scmiclassical face, this is in fact the same construction as in [IB], but with the 
parameter h. Using equation ((5 . 3|) and the projective coordinates (|5.5p . we find 
that the operator P(h, a) lifts to 

„2 

P (h, a) = h 2 {-{Xd x ) 2 + nXd x + P AX 2 8 X + X 2 A Hc + p 2 X 2 W - — ) - a 2 . 

In these coordinates the left face is given by {X — 0}. Therefore, the kernel of the 
composition K(P(h, a)G) when lifted to B n+1 Xq B n+1 is, near the left face, equal 
to 

K(P(h, o)G) = (h 2 {-{Xd x ) 2 + nXd x ~\)~ ^) K (G) + 0(X 2 ). 

The operator I{P{h,<r)) = h 2 (-(Xd x ) 2 + nXd x - x) - is called the indicial 
operator oiP(h,a). Since G G $~^ +i ^ , -a _1 '^ +, '*(B n + 1 ), then near the left face 
K{G) G ?C2+if -f-b?+if (Mq, r ). But I(P(h,a))X% +i % = 0. So we deduce that 
near C, K{P{h,a)G) G K.% +i % +1 >~% (M ,r). That is, we gain one order of 
vanishing at the left face. Since we already know from (|5.14j) that the kernel of the 
error vanishes to infinite order A and at the front face J 7 , and x — Rpn, x 1 — Rpl, 
the kernel of the error R, on the manifold B™ +1 x B n+1 satisfies 

K{R) G /^x^+^ + V t+if C°°(B" +1 x B™ +1 ). 

Then one can use a power series argument to find G4 with Schwartz kernel in 
h 00 x^ +i i +1 x'^ +l ^C 00 (M n+1 x B" +1 ) such that 

P(h, <j)G 4 -Re $-~.~.oo,t-Mf (B" +1 ). 
So G = Go - Gi + G 2 - G 3 - G 4 G . g thg desired parame trix. 

6. L 2 -BOUNDS FOR THE SEMICLASSICAL RESOLVENT 

We now prove bounds for the semiclassical resolvent, R(h,a) = P(h,cr)~ l : 

Theorem 6.1. Let M > 0, h > and a G C be such that ^f- < M. Let a,b > 
max{0, ^ff~}- Then there exists ho > and G > independent of h such that for 
h G (0,h ), 

(6.1) \\x a R(h,a)x b f\\ mKn+1) < Ch-'-^WfW^^y 
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As usual, we prove Theorem 16.11 by obtaining bounds for the parametrix G 
and its error E on weighted L 2 spaces. As a preliminary remark, we recall that 
elements of VP2 on compact manifolds without boundary are £ 2 -bounded with an 
/i-independent bound; the same holds for elements of h . Thus, the diagonal 
singularity can always be ignored (though in our setting, due to negative orders of 
the operators we are interested in, Schur's lemma gives this directly in any case). 

The following lemma follows from the argument of Mazzeo [HI Proof of Theo- 
rem 3.25], since for the L 2 bounds, the proof in that paper only utilizes estimates 
on the Schwartz kernel, rather than its derivatives. Alternatively, it can be proved 
using Schur's lemma if one writes the Schwartz kernel relative to a b-density. 

Lemma 6.2. Suppose that the Schwartz kernel of B (trivialized by \dgs(z')\) sat- 
isfies 

\B{z,z')\<Cplp P R , 

then we have four situations: 

If a, /3 > n/2, then \\B\\c(v>) < C'C. 
Ifa = n/2, /3>n/2, then \\\logx\~ N B\\c(l 2 ) < C'C, for N > ^. 
If a > n/2, P=n/2, then \\\B\ log x\- n \\ C (l') < C'C, for N > 1 
If a = (3 = 71/2, then\\\\ogx\- N B\\ogx\- N \\ c(L 2 ) <C'C, N > ^ 

Now let B{h,a) have Schwartz kernel B(z, z' ,o~,h) supported in r > 1, and 
suppose that 

B(z,z',a,h) = e-™ r / h h k p n L /2+ ' ( p n B /2 x a (x') b H, H g L°°, and ^ < N. 

Since from Proposition 12.71 r = —logpRpL + F, F > 0, e~ t ^ r = p % ^ e~ l ^ F , 
and ^ <N, \e-^ F \ < C = C(N), it follows that 

\B(z,z',a,h)\ < Ch k pl /2+1 - llaa/h+a p n R /2 - Ila<T/h+b R a + b 

As an immediate consequence of Lemma l6.2[ if a+b > 0, 5q > and 7— lma/h+a > 
Sq and — Imcr /h + b > 5q, then 

\\B\\ L 2 < C'Ch k . 

If either 7 — lma/h + a = or — Jmcr/h + b = 0, we have to add the weight 
I \ogx\~ N , with N > 5. On the other hand, suppose now that the Schwartz kernel 
of B is supported in r < 2, and (again, trivialized by \dgs(z')\) 

\B(z,z',a,h)\ < Ch k (r/h)- e (h/r) s , s < n + 1. 

Note that for fixed z',a,h, B is L 1 in z, and similarly with z' and z interchanged. 
In fact, since the volume form is bounded by Cr n dr{z') du)' in r < 2, uniformly in 
z, Schur's lemma yields 

f 2 

\\B\\ L 2 < CC"h k / (h/r) s {r/h)- e r n dr. 
Jo 
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But 

(h/r) s (r/hy l r n dr 

< C (h/r) s r n dr + J (r/h)- e r n di^j = d{h n+1 + h' e ), 
so we deduce that 

\\B\\ L 2 < CC'{h n+1+k + h k " e ). 
First, if, with n + 1 = 3, G is the parametrix, with error E, constructed in 
Section 21 then, writing G = G\ + G 2 , and provided ^jf- < 1 and \a\ > 1 (recall 
that Ui has a factor with G\ supported in r < 2, G2 supported in r > 1, 

then I Gi| < Ch- 2 r-\ \G 2 \ < Ce lnl,yr / h p^ 2 ' p n J 2 , and thus, using Proposition [23 

|:r Q (2/) b Gi(z,z',<x,/i)| < Ch^- 1 (h/r)^ 1 (r/h)- 11 ' 2 , 
\x a (x') b G 2 (z,z',a,h)\ < Ch- n - 1 x a (x , ) b p n L /2 - lm ' 7/h p n R /2 - lm ' T/h . 
On the other hand, 

(6.2) \x a (x') b E(z, z', a, h)\ < Chx a (x') b p n L /2+2 - Ima/h p n J 2 - lma/h . 

Thus, we deduce the following bounds: 

Proposition 6.3. Suppose n + 1 = 3. Let G(h,o~) be the operator whose kernel is 
given by (gJ2), and let E{h, a) = P{h, a)G(h, a) - Id. Then for \a\ > 1, ^f- < 1, 
a > a > 0, and ^f- < b < 2 - b > 0, we have, with C independent of h, 

\x a G(h,a)x b 7||l2(b»+i) < C/i -1- ^ I |/|| i2(B „ +1) an d 

\\x~ b E(h, a)x b f\\ L 2 (Kn+ i ) < Ch\\f\\ L 2(nn+iy 

If either a — ^jf- or b — ^jf- , or a = b = ^jf-, one has to replace the factor x^~ ~ 

, Imt7 \±i 
in (|6.3|) with I x"^ 2 "! \ogx\~ N j , iV > i, to obtain the L 2 bounds. 

In view of (|6.2[) this cannot be improved using the methods of section [4] To 
obtain bounds on any strip we need the sharper bounds on the error term given 
by Theorem 15.11 We now turn to arbitrary n and use the parametrix G and error 
E, constructed in Theorem OD] Writing G = G + Gi + G 2 , with G G G x 
supported in r < 2, G 2 supported in r > 1, Gi € e-W'ij/"^ ' 00 ' a- -1 ' 00 ^ g< 2 g 
g-iar/^-^.f +2 then for | CT | > 1 and ^ < JV, 

x a G Q (z,z',a, h){x') b e%% 

\^{^) h aL{z,z\a,h)\ < Ch^ 1 {r/h)- n / 2 , 

\x a {x') b G 2 (z 1 z' 1 a,h)\ < Ch- n - 1 x a (x') b p n L /2 - IlrL ' y/h p n J 2+Inia/h . 

On the other hand, writing E = E\ + E 2 , with E\ supported in r < 2, E 2 supported 
in r > 1, Ei G h^p™^™^ ~ 1,(X> , E 2 G h°° pf^^^^^ , then for any fc 
and M (with G = G(M, TV)), 

Ix- 6 ^') 6 ^!^,/,^/!)! < Ch k , 

\x- b (x') b E 2 (z,z',a,h)\ < Ch k x M (x') b p n B /2 - Ima/h . 
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In this case, we deduce the following bounds: 

Proposition 6.4. Let G(h.a) be the operator whose kernel is given by (|5.ip . and 

< m, M > 0, a > ^ 



let E(h,cr) = P(h,a)G{h,a) - Id. Then for \a\ > 1, ^ < M, M > 0, a > Hs., 



a > 0, and b > ^jf~-, b>0, and N arbitrary, we have, with C independent of h, 

\\x a G(h,a)x b f\\ L 2 (Kn+ i ) <Ch- 1 -^\\f\\ L 2 {Bn+ i ) and 
"'' 4 ' \\x- b E(h,a)x b f\\ LKMn+ t } < Ch N \\f Wvpn+iy 



If either a = ^jf- or b = ^r 2 -, or a = b = ^- 2 -, one has to replace the factor x :iz ~ B ^ L 

/ ImtT 1 \±i 

in (|6.4|) with ( x~ ~ | log x\~ k ) , A; > \, to obtain the L 2 bounds. 

Now we can apply these estimates to prove Theorem 16.11 We know that 
P(h,a)G(h,a) = I + E(h,a). 
Since R(h, a) is bounded on L 2 (B n+1 ) for Imcr < we can write for Im a < 0, 

G{h, a) = R(h, a)(I + E(h, a)). 
Therefore we have, still for Imcr < 0, 

x a G(h, a)x b = x a R{h, a)x b {I + X - b E{h, o-)x~ b ). 
For a, b and a as in Proposition 16.41 we can pick ho so that 

\\x- b E{h,a)x h f\\ L ^ L , < \. 

In this case we have 

x a G(h, a)x b {I + x- b E(h, a) X - b )- 1 = x a R(h, a)x b 
and the result is proved. 

7. Proof of Theorem 11.11 

Now we are ready to prove Theorem 11.11 To avoid using the same notation 
for different parameters, we will denote the spectral parameter in the statement of 
Theorem 1 1.1 1 by A, instead of a. We write, for | Re A| > 1, 



1 '* + sw*-\- * 



(Re A) 2 V 4 J (Re A) 



Thus, if we denote h = and a — -^j, and E(h, a) and G(h, a) are the operators 
in Proposition ^. 41 we have 

< a » + ^ - * - T^ity isa> = ( M +£ <5^ ikx 

Since i?a(A) = (^A gs + x 2 W — A 2 — is a well defined bounded operator if 

ImA < 0, we can write, 

G ^R^) = ( RcA ^^( Id+ ^'^)' ImA<0 ' 
Therefore, 
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I^tJTT'^^/IIl^+i) < C(ReA) 1+ t||/|| L2(B „ +1) . 



According to Proposition 16.41 if Im A < M, we can pick K such that if | Re A| > K, 
and Im A < b then, 

''^R^R^''4 

Therefore (Id +x~ b E( -^rj, rZx)^) 1 ' 1S holomorphic in ImA < M, and bounded 
as an operator in L 2 (B 3 , g), with norm independent of Re A, provided b > ImA. On 
the other hand, if a > ImA, then from Proposition! 

1 A_ 

v ReA' ReA ; 

Since, 

x a Rs(X)x b = (Re rV^, (/ + X -»E(^- y ^ ' 

then, for and for a, b in this range, and | Re A| > K, x a R s (X)x b is holomorphic and 

\\x a R s (X)x b f\\ L 2 (v ^ g) < C(Re \)*- 1 \\f\\&<p>, a y 

When either a — Im A or b = Im A we have to introduce the logarithmic weight 
and in Proposition 16.41 This concludes the proof of the L 2 estimates of Theorem 

o 

The Sobolev estimates follow from these L 2 estimates and interpolation. First we 
observe that the following commutator properties hold: There are C°° (B 3 ) functions 
Ai and Bj, i = 1,2, and 1 < j < 5, such that 

[A gs ,x a ] = A 1 x a +A 2 x a xD x , 

[A g5 ,x a (\ogx)- N ] = B lX a (logx)- N + B 2 x a (logx)- N - 1 + 

B 3 x a (logxy N ~ 2 + (BiX a {logx)- N + Bsx^logx)-"- 1 ) xD x . 

Hence 

A gs x a R s (a)x b v = x a A gs R s {a)x b v + A lX a R & {a)x b v + A 2 x a xD x R s {a)x b v. 

Since a. b > 0, A gs is elliptic, and A gs Rs(cr) = Id + (cr 2 + 1 — x 2 W(x))R$ it follows 
that, see for example [IS], that there exists a constant C > such that 

\\x a Rs(a-)x b v\\ H 2^ 

(7 1) 

< C ^cr 2 ||x b u|| L 2( B 3) + \\x a Rs(cr)x b v\\ L 2 {B 3 ) + \\x a R s (a)x b v\\ H i {B 3- ) 
By interpolation between Sobolev spaces we know that there exists C > such that 
(7-2) IK«IIh i ( b3) < C\\x a v\\ L 2 {B3) \\x a v\\ H 2 m 

Therefore, for any e > 0, 

(7-3) \\x a v\\ H i (B3) < C (e||a; a w|| ff 2 (B3) + e _1 ||x°i;||£2(B3)J , 

and if one takes e small enough, (|T. 1 [) and (|T.3[) give (| 1 .3() with k = 2. If one uses 
(|TT2|) and (TQ|) with fc = 2 one obtains (TO)) with k = 1. The proof of (fl~5j) follows 
by the same argument. 

This completes the proof of Theorem ll.il 
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8. Structure of Ax near the boundaries 
We begin the proof of Theorem |4] by analyzing the structure of the operator Ax 

2 dr l 



near r = rbH and r = r s i. We recall that /3(r) = \^ct 2 {r) and that /?bH = /^mi) 



and p s i = /3(r s i). 

We show that near these ends, after rescaling a, the operator a% Axa~^ is a 
small perturbation of the Laplacian of the hyperbolic metric of constant negative 
sectional curvature — /3 2 H near rbH and — near r s \. 

Since a'(r) ^ near r = rbH and r — r s i, Ax can be written in terms of a as a 
'radial' coordinate 

A x = f3r~ n aD a (pr n aD a ) + a 2 r" 2 A w . 

We define a C°° function x on [rbH,?"si] which is positive in the interior of the 
interval and rescales a near the ends by 

(8.1) a = 2rbH~A>H2; near r = rbH and a = 2r s i|/3 s i|x near r — r s i- 

Using x instead of a near the ends rbH and r s i, we obtain 

A x = Pr' n xD x {pr n xD x ) + 4x 2 ^r^r- 2 A u , near r = r bH , 
A x = /3r- n xD x (pr n xD x ) + Ax 2 ^r^r' 2 A Ul near r = r sI . 



Proposition 8.1. There exists S > such that, if we identify each of the neighbor- 
hoods of {x — 0} given by r 6 [rbH! r bH+5) andr G (r s j — (5, r s j], toztft. a neighborhood 
of the boundary of the ballW l+1 , £/ien i/iere ea;is£ two C°° functions, Wbii(x) defined 
near r = rbH, and W s i(x) defined near r = r s i, such that 

a?Axa~2 — x~ Axx~~ — A gm + x 2 WbH — 0bH~T> near r — r bH and 
(8-3) > \ 2 

cfi AxaT^ — x% Axx~% = A gsI + x 2 W s i ~ Pli—, near r = r s i, 

where gbH and g s i are small perturbations of the hyperbolic metrics with sectional 
curvature ~Pl H and —/3 2 j respectively on the interior o/B n+1 , i.e. 

Adz 2 Adz 2 

(8 - 4) 9bH = p 2 bH {i - M 2 ) 2 + HhH ' and 9sI = fU^VA 2 ) 2 + Hsh 

where HbH and H s j are symmetric 2-tensors C°° up to the boundary o/B n+1 . 



Proof. It is only necessary to prove the result near one of the ends. The computation 
near the other end is identical and one only needs to replace the index bH by si. 
From (|8.2j) we we find that near r = rbH, 

zi A x x-2 =/3 2 (xD x ) 2 + m/3 2 xD x + f3 r - n (xD x ((3r n ))xD x + {2/3 hH r h}i ) 2 r- 2 x 2 A w 

2 

-z|/3r-"^(/3r")-/3 2 ^. 
Let <?bH be the metric defined on a neighborhood of dM n+1 given by 

(8-5) ^h = ^2^2 + A bH r2 "^2~' where A b H = 2|/? b H|rbH- 

The Laplacian of this metric is 

A 9bH = P 2 {xD x ) 2 + in/3 2 xD x + /3r- n (xD x (/3r n ))xD x + X 2 K r- 2 x 2 A UJ . 
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Therefore we conclude that near the ends r = rbH 

2 

x% A x x~% = A 3bH - f3 2 ^ - i^f3r- n xD x (f3r 2 ). 

Since r = r(x 2 ), we can write near and r s i, 

r = r hu + x 2 A hii (x 2 ) and /3(r) = /3 bH + x 2 B hH (x 2 ) near r = r hli . 
Therefore, near r = 

= -^+x 2 B hH (x 2 ). 

We conclude that there exist a symmetric 2-tensor H\xi(x 2 , dx, duj) near r = rbH 
which is C°° up to {a; = 0}, and such that the metric c/bH given by (|8.5[) can be 
written near r — rbH as 

dx 2 doj 2 

(8.6) 5bH = -S2—2 + AR 2 r 2 + H ™ near r = r bH- 

PbH x *PbH x 



Let g be the metric on the interior of B n+1 which is given by 

. _ A\dz\ 2 
9 ~ c 2 (l - \z\ 2 ) 2 ' 

We consider local coordinates valid for \z\ > given by (x,uj), where uj — z/\z\ 
and x = 1+ L ■ The metric g written in terms of these coordinates is given by 

dx 2 . 9n9 duJ 1 
9 = 3-2 + (1 - x 2 ) 2 



4c 2 a; 2 

Therefore, near x = 

dx 2 duj 2 2 , , . 

.9 = -o-o + 7TT + ^(^ ,u,dx,dw), 
c x Ac x 

where H is a symmetric 2-tensor smooth up to the boundary of B n+ . This concludes 
the proof of the Proposition. □ 

9. From cut-off and models to stationary resolvent 

Next we use the method of Bruneau and Petkov [3] to decompose the operator 
R(<t) in terms of its cut-off part x-R(c)x and the contributions from the ends, which 
are controlled by Theorem ll.il For that one needs to define some suitable cut-off 
functions. For 5 > let \ji Xji and Xji 3 = 1, 2, defined by 



xiO) = 


= 1 if r > r bH + 


4<5, 


Xi(«) = 


if r < - 


f 35, 


X\(r) 


= 1 if r > r b H ^ 


-25, 


Xl(a) = 


= if r < r bH 


+ <5, 


Xi(r) -- 


= 1 if r > r bH + 


66, 


Xi(a) = 


if r < r b H - 


f 5(5, 


X2O) 


= 1 if r < r s i — 


46, 


X2{a) = 


if r > r s i — 


3(5, 


X\(r) 


= 1 if r < r s i - 


-26, 


x\{oi) = 


- if r > r s i - 


-5, 


Xi{r) 


= 1 if r < r s i — 


66, 


X2{a) = 


if r > r s i - 


5(5, 



and let 



Xs(r) = 1 - (1 - Xi)(l - Xi) - (1 - X 2 )(l - X2)- 
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X3(r) is supported in [rt,H + S, r s i — S] and Xi{ r ) = 1 if r € [rbH + 26, r s i - 25}. Let 
X G C5 (r bH , r s i) with X (r) = 1 if r £ [r bH + 5/2, r a ~ 6/2]. 

Now we will use Proposition 18.11 for 5 small enough. Let <7bH and g s i be the 
metrics given on the interior of B™ +1 given by (|8.4[) and let 

(9.1) 

Adz 2 Adz 2 
9bH,8 = ^ _ | z | 2 )2 + ( X - Xi)^bH, and 3sU = _ + (1 - x 2 )H sl . 

Since gbn,s = ffbH if Xi = 0> an d g s i,s = <7si if X2 = 0, it follows from Proposition 
18.11 that, for x given by equation (|8.1|) 



a^A x a-T(i _ = (A 9bH _ 4 + x 2 VK bH - — /?bH - <7 2 )(1 - Xi)/, and 

o 

^Axa^(l - Xi)/ = (A 3bH . a + .x 2 VK bH - — /? b H - <7 2 )(1 - Xi)/- 

Q 

« T A**" t (1 - X2 )f = (A 9sM + z 2 W sI - T /? s 2 i - a 2 )(l - X2)L and 
at Axa't (1 - x i)/ = (A 3sM + x 2 W sI - ^/3 2 - <7 2 )(1 - X J)/. 

Let 



= a 2 R{a)a 2 =(a 2 A x a 2 - a 2 ) , Im cr < 0, 

and let 

n 2 

R hH (a) = (A 5bH , + x 2 V7 b H - o 2 - T /3 2 H ) 1 and 

2 

i? sI (a) = (A SsM + x 2 W sl -a 2 - -fey 1 

be operators acting on functions defined on B' i+1 . 

If, as in Proposition 18. 11 we identify neighborhoods of r = rbH and r = r 8 i with 
a neighborhood of the boundary of B™ +1 , we obtain the following identity for the 
resolvent 

R a (a) = R a (a)x 3 + {l-Xi)Rbu((r)(l-xl) + ^-X2)R s i((T)(l-xl)- 

(9.2) - 1 

R a (a)lA^ u ,l-xi]RbH(<r)(l-Xi)-M<r)[\ei^-X2}R s i(<T)(l-xl)- 
Similarly, one obtains 

R a (a) = xM*) + (1 - Xi)-RbH(a)(l - X i) + (1 - X2)^iW(l - X 2 )+ 

(9.3) , 1 
(l-xl)Rhn((r)[A gb3 ,l-xi]R a ((T) + (l-xi)R s i^)[^A-X2]Ra^)- 
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These equations can be verified by applying a? Ax a 2 — a 1 on the left and on the 
right of both sides of the identities. Substituting (|9.3I) into (|9.2[) we obtain 

R a (a) = M 1 (a)xR a (<j)xM 2 (<j)+ 

(1 - Xi)Rbn(<r)(l - x \) + (1 - X 2 )i?si(a)(l - xl), 

where 

(9.4) Af 1 (o-)=X3 + (l-Xi)-RbH(o-)(l-Xi)[A ffbH ,l-Xi] 

+ (l- x i)ii sI (a)(l-X2)[A gsI) l-X2], 
M 2 (a) = X 3 - [A 9bH , 1 - xi](l - Xi)JM<r)(l ~ Xi) 

-[A Ssr ,l-X2](l-X 2 Ri(0(l-X2). 

This gives a decomposition of R a (a) in terms of the cutoff resolvent, studied 
by Bony and Hafner PQ , and the resolvents of he Laplacian of a metric which are 
small perturbations of the Poincare metric in B n+1 . The mapping properties of such 
operators were established in Section [7] Next we will put the estimates together 
and finish the proof of our main result. 

10. The proof of Theorem 0] in 3 dimensions 

We now prove Theorem |4] for n + 1 = 3 using Theorem 11.11 and Theorem [3] We 
first restate a strengthened version of the theorem which includes the case where 
the weight b = Im a. 

Theorem 10.1. Let e > be such that (fit))) holds and suppose 

0<7<min(e,/3 fcff , |/3 S/ |,1). 
Then for b > 7 there exist C and M such that iflraa < 7 and | Re<r| > 1, 

(10.1) \\a b R{a)a b f\\ L2{x , n) < C|(7| M ||/|| i2(X;n) , 

where a was defined in (|15p and the measure f2 was defined in (JSj) . Moreover, for 
N > 5 and < S « 1, choose ipN( r ) £ C°°( r &.£r, r s j) with ipjyfr) > 1, such that 

(10.2) ipN = I loga| _iV «/ r — rbH < <5 or r s j - r < (5. 

TTien wra£/i 7 as above, there exists C > and M > smc/i that for | Reu > 1 and 
Im cr I < 7 

||a Im ^(a)i?(a)5 b /|| i2(X;n) < C|a| M ||/|| LW) , 

(10.3) \\a b R(a)a Im ^ N (a)f\\ L 2 (Xin) < C\a\ M \\f\\ L 2 (Xin) and 
\\a lm ^ N (a)R(a)a lm ^ N (a)f\\ LHX , n) < C\a\ M \\f\\ LHX .^. 

Proof Recall from (l9Tl) that 

1 , 1 

gbn,5 = -fi2~9s and g s i,s = ^r.9<5, 
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where gs is of the form (|1.2p . So we obtain A 9bH s — (3^ H A g5 and similarly A 9aI 4 



P%A gs . Therefore, 



(10.4) 



n 2 

Rbn(v) = ( /3b H A ss + x 2 W hu -a 2 - —p 2 H 



(3-4 A gs + x'P^W hn - a 2 ^ - - = fi^RM^V 1 ) 



-l 



Therefore, by replacing a with cr | /?bH I ~ 1 in (|l-3[l and (|1.5I) setting a — p^j, and 
6 = |^ h | , we deduce from Theorem 11.11 that there exists Sq > such that if 
< S < S , for Imcr < A and Im a < B and | Re tx| > K(S), 

R hH (a)xT^ v\\ H k {m3) < C\a\ k \\v\\ L 2 (B3)l k = 0, 1,2, 
(10-5) _g_ 

Hll^HlflbHCT)! lObHuH^BS) < C|f7| fe ||w|| H -, (B3) , fc = 0,1,2. 

Of course, the same argument applied to R s i gives 

| x ~\ Rsi (a)x Tfcfl v 1 1 H k ( B 3 j < C|a| fc || V || L 2 (B3) , fc = 0,1,2, 

(10.6) _g_ 

||a:^ni? sI( r)a;^ru|| L 2 (B3) < CVplM^-*^, fc = 0,l,2. 

When A = Imcr, or B = Ima, we define Thn,A,B,N and TLi,A,b,jv an d as in (|1-4|) 
by replacing Rs(cr) with either i?bH(c) or J? s i((t), a with ^4 and & with B. Using 
(|1.5[) we obtain for J = bH or J = si, 

\\Tj,A,B,N(v)v\\ H k (B3) < C|CT| fe ||t;|| L 2 (B 3 ) , k = 0,1,2, 

(10.7) 

||Tj,A,s,iv«||z, 2 (B3) < C|(t| ||u|| h -»(b3), fe = 0,l,2. 

Now we recall that a = 2rbH/3bH near rbH and similarly a — 2r s i/3 s i near r s i. We 
will use these estimates, identity (|9.4j) and Theorem|3]to prove Theorem[4l Indeed, 
in the case a > Im a, b > Im a we write 

a a R a a b = a a M 1 {a)a b a- b xR a {<j)xa- b a b M 2 (cr)a b + 

(1 - X i)5 a i?bH(cr)a b (l - Xl ) + (1 - X2 )a a i? sI (cr)a b (l - X \)- 

Notice that the measure in B™ +1 is x~ n dxdu>, which corresponds to a~ n ~ 1 dadui 
which in turn corresponds to a~ n ~ 2 drduj. In this case n = 2, but this part of the 
argument is the same for all dimensions, and we will not set n = 2. Thus, we deduce 
from Theorem 11.11 that 

11(1 - Xl)® a Rbli(cr)a b (1 - Xl)\\L 2 (X;a-"-- 2 drdui) < C\ I V 1 1 L 2 ( X;a- n - 2 drdu) > 
||(1 - X2)a a J Rsl(cr)a h (l - X^HU 3 ^-^-^^) < C| |u| | i 2 (X;Q -„-2 (ir(iw ). 

Recall that i7 = a~ 2 r 2 drdu). Since r £ [rbH, r s i], 7"bH > 0, this gives 

IK 1 - Xi)a a &~^ Rbn((r)a% & b {l -xi)IU=(x ; n) < C\\v\\ L i {x . a) , 

(10.8) _„ „ 

\\(l- X2 )a a a 2 R s i{cr)a 2 - X ^)v\\ L 2 (x . n) < C\\v\\ L 2 (x . a -n-2 n) . 

Similarly, using the Sobolev estimates in Theorem 1 1.1[ we obtain 



a a a 2 Mi{a)a 2 a b v\\ L 2 [x . n] < C\a\\\v\\ L 2 {x . n) , 
a a a~% M 2 (cr)a^ a b v\\ L 2 {x . n) < C\a\\\v\\ L 2 {x . n) . 
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Since \ is compactly supported in the interior of X, it follows from Theorem [3] 
that 

(10.10) \\ar b xRa{o)x^ b v\\L^(X:n) < C\\v\\ L 2( X , n) . 

Estimates (fTM]) . (fT09l) and (|10.10|) imply that 

\\a^ b a^^ R a {(j)a^ x^~ bv \\L^(x-n) < C\\v\\ L 2 (x . n) . 
But a-%R a (a)a% = R(a). This proves (fTTJTTj) . □ 

11. The proof of Theorem H] in general dimension 

We will outline the main steps necessary to connect to the results of [5], and 
refer the reader to [5] for more details. 

First choose S so small (IT21) holds. Let Xo and Xi be as defined in (fTTj). then we 
recall that for <5 small, 

or^Pa~^\ Xo = P(h,a), 

where P(h,a) stands for model near either rt,H or near r s \. By Theorem 16.11 there 
exists ho > such that for h 6 (0, /lo), 

||a; Q (fo 2 P -o-) -1 ^ 6 !!^^ < C/r 1 "* ere (1 - c, 1 + c) X (-c,eh) C C. 

Let Pi be the operator defined in (IT5| and let e > be such that (TH|) holds. 
Then it follows from Theorem 2.1 of |S] that there exist hi > 0, C > and X > 
such that for h G (0, /zi), 

||5 h a-*(/i 2 A x -cr 2 )- 1 a Q a _ t|| i2 L2 < C7i -K , 
(11.1) — 

(7€(l-c,l + c)x(-c,e/i)cC. 

The estimate in Theorem |4] follows by restating (|11.1|) in the non-semiclassical 
language, i.e. multiplying it by h 2 and replacing a by h~ 1 a. 
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